1 L ..  P[  rCi*<P  4.  I ItON 

,OAv#a  vzhool 

MoMbnat  CAUK^N)^ 

(\oo<*  Sts 


RMC 

Research 

Corporation 

m A Resource  Management  Corporation  Subsidiary 


The  RMC  Systems  Group 


Report  UR-22G 


IMPLEMENTATION  OF  A MODEL  OF  REQUISITION  PROCESSING 
FOR  THE  SHIPS  SUPPLY  SIP* PORT  STUDY 


Project  Director:  George  M.  Lady 
Cbnsultant:  Carl  M.  Harris 


October  23,  1973 


BEST 

SCAN 

AVAILABLE 


Prepared  for 

Naval  Supply  System  Command 


7910  Woodmont  Avenue.  Bethesda,  Maryland 20014 


RMC  Report  UR-226 


IMPLEMENTATION  OF  A MODEL  OF  REQUISITION  PROCESSING 
FOR  THE  SHIPS  SUPPLY  SUPPORT  STUDY 


Project  Director : George  M.  Lady 
Consultant:  Carl  M.  Harris 


October  23,  1973 


Prepared  for 

Naval  Supply  System  Command 


sy-% 


CONTENTS 


1.  INTRODUCTION 1-1 

Modeling  Requisition  Processing  1-2 

Capacity  Factors 1-4 

Outline  of  This  Report 1-5 

2.  METHODOLOGY  2-1 

Steps  in  the  Methodology 2-3 

Step  I:  Eliminate  Cycles  from  Flow  Matrix  2-5 

Step  II:  Establish  Work  Load  Factors 2-5 

Step  III:  Establish  Capacity  Factors 2-6 

Step  IV:  Implement  Multipriority  Queueing  Model 2-10 

Step  V:  Path  Analysis 2-10 

Step  VI:  Establish  Frequency  Distribution  of  Throughput 

Times  Per  Path 2-11 

Step  VII:  Calculate  Distribution  of  System  Throughput  Times  . . 2-11 

Illustration  of  the  Methodology  2-12 

Step  I:  Eliminate  Cycles 2-12 

Step  II:  Establish  Work  Load  Factors 2-13 

Step  HI:  Establish  Capacities  2-15 

Step  IV:  Implement  Multipriority  Queueing  Model  2-16 

Steps  V and  VI:  Path  Analysis  and  Path  Frequency  Distributions.  2-23 

Step  VII:  System  Throughput  Times 2-29 

3.  MODELING  REQUISITION  PROCESSING  AT  SPCC  3-1 

Step  I:  Eliminate  Cycles  from  Flow  Matrix  3-3 

Step  V:  Path  Analysis 3-4 

Step  III:  Establish  Capacity  Factors 3-7 

Step  II:  Establish  Work  Load  Factors 3-9 

Step  IV:  Implement  Multipriority  Queueing  Model 3-9 

li 

♦ cat 


Outputs  of  the  Multipriority  Queueing  Model  3-16 

Step  VI:  Establish  Frequency  Distributions  of  Throughput 

Times  Per  Path  3-25 

Path  Enumeration 3-31 

4.  CONCLUDING  REMARKS 4-1 

APPENDIX  A A-l 

APPENDIX  B - Program  Listings B-l 


1 


INTRODUCTION 

This  report  presents  the  results  of  implementing  a methodology  to  relate  Navy 
Supply  System  throughput  time  associated  with  requisition  processing  and  ma- 
terials handling  to  the  resources  dedicated  to  those  activities.  This  effort  on 
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the  part  of  the  RMC  Research  Corporation  is  in  support  of  the  S (Ships  Supply 
Support  Study)  project  which  is  a general  simulation  of  supply  system  activities. 
Consideration  of  alternative  methodologies  and  selection  of  a particular  method 
were  accomplished  by  RMC  Research  Corporation  personnel  prior  to  the  study 
effort  reported  on  in  this  document.  * 

In  developing  the  methodology,  supply  system  activities  were  studied  in 
the  areas  of  communications,  transportation.  Inventory  management,  and  stock 
point  operations.  Pilot  applications  of  the  method  and  other  data  revealed  the 
throughput  times  associated  with  communications  activities  were  not  large  enough 
to  warrant  an  expensive  modeling  effort.  Transportation  activities  could  be  suf- 
ficiently modeled  on  the  basis  of  an  enumeration  of  transport  vehicles,  schedules, 
and  various  administratively  determined  policies.  The  interesting  and  seemingly 
important  portions  of  throughput  time  for  modeling  purposes  were  associated  with 
requisition  processing  and  materials  handling  at  inventory  control  points  and  stock 
point  facilities. 

Further  study  of  materials  handling  operations  revealed  that  resource  level/ 
throughput  time  relationships  could  be  reasonably  explained  by  simple  queueing 


1.  See  RMC,  Inc.  Final  Report,  UR-176,  Methodology  for  the  Measurement 
of  the  Relationship  Between  Naval  Supply  System  Resources  and  Supply  System 
Throughput  Time,  June  5,  1972. 
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theoretic  formalisms  and  an  incorporation  of  pickup  and  delivery  schedules. 
Generally,  the  data  describing  materials  handling  operations  revealed  very 
modest  backlogs.  Confidence  in  the  throughput  times  predicted  by  an  extrapo- 
lation of  materials  handling  models  was  enhanced  by  the  fact  that  the  models 
could  be  designed  to  closely  replicate  the  materials  handling  process.  Only 
that  portion  of  throughput  time  associated  with  requisition  processing  appeared 
to  offer  complexities  which  would  not  submit  to  main  force  analytic  techniques. 

MODELING  REQUISITION  PROCESSING 

Study  requisition  processing  at  inventory  control  points  (ICPs)  and  stock 
points  revealed  long  waiting  times  not  apparently  related  to  servicing  the  requi- 
sitions. Certain  of  these  times  could  be  explained  by  policy  and  administrative 
procedures  not  related  to  resource  levels.  However,  in  the  initial  modeling 
efforts  waiting  times  due  to  backlogs  could  not  be  explained  by  the  relatively 
simple  queueing  theory  models  attempted. 1 Several  factors  could  explain  the 
failure  to  model  backlogs: 

• the  queueing  models  attempted  did  not  take  into  account  the  system  of 
priorities  under  which  requisitions  were  processed, 

• the  resources  dedicated  to  requisition  processing  are  often  used  in  a 
number  of  other  capacities, 

• the  assumptions  about  work  load  arrival  rate  and  service  time  distribu- 
tions used  by  the  queueing  model  might  be  unrealistic,  and 

• the  backlogs  might  be  contrived  by  supply  system  personnel  for  any  of 
a number  of  reasons  ranging  from  a desire  to  efficiently  use  personnel 
to  foster  job  security. 

The  emphasis  of  the  methodology  implementation  reported  on  here  was  to 
more  intensively  study  and  model  requisition  processing  especially  at  inventory 
management  facilities.  A serious  difficulty  was  the  lack  of  a suitable  supporting 
data  base.  For  reasons  of  economy  the  study  was  constrained  to  secondary  data 
sources.  An  investigation  of  information  contained  in  the  Military  Standard 


1.  A highly  detailed  simulation,  SIMCOM,  conducted  by  FMSO  personnel  to 
model  requisition  processing  at  SPCC  also  failed  to  predict  waiting  times  due  to 
backlogs . 
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Requisitions  and  Issue  Procedures  (MIIiJTRIP)  reporting  system  revealed  that 
sufficient  data  were  available  although  such  data  were  not  currently  reported  in 
a usable  fashion.  The  automatic  data  processing  arrangements  at  both  Ships 
Parts  Control  Center  (SPCC)  and  Electronics  Supply  Office  (ESO)  maintain  what 
is  called  a ’’requisition  status  file."  The  purpose  of  this  file  is  to  enable  the  in- 
ventory management  facility  to  determine  the  location  (in  its  administrative  struc- 
ture) and  status  of  a requisition  as  needed,  usually  in  response  to  an  inquiry  from 
the  originator  of  the  requisition.  Included  among  the  data  contained  in  this  file 
are  the  "date  of  receipt"  by  the  facility,  the  particular  branch  or  division  of  the 
facility  in  which  the  requisition  is  located  at  a given  point  in  time  (the  "local 
routing  code"),  the  "date  of  last  action"  which  usually  denotes  its  initial  receipt 
in  its  current  location,  and  its  status  (e.g. , under  manual  review,  on  back  order, 
under  procurement,  etc. ).  The  sum  total  of  these  data  would  document  the  fol- 
lowing: 

• the  particular  administrative  subdivisions  of  the  facility  which  serv- 
iced a given  requisition, 

• the  numbers  of  requisitions  received  by  each  subdivision  per  day, 

• the  distribution  of  times  a requisition  spent  in  each  administrative 
subdivision, 

• the  frequency  with  which  requisitions  experienced  alternative  final 
actions  (e.g. , the  percent  sent  to  stock  points,  placed  on  backorder, 
placed  under  procurement,  etc.,  . . .), 

• the  frequency  with  which  requisitions  received  by  a given  subdivision 
were  sent  on  to  each  other  subdivision,  and 

• the  distribution  of  total  times  from  initial  receipt  of  the  requisition 
until  its  final  disposition. 

The  problem  with  the  status  file  records  as  kept  was  that  the  prior  "history"  of  a 
requisition  was  erased  when  the  file  was  updated. 

An  investigation  by  Fleet  Material  Support  Office  (FMSO)  personnel  revealed 
that  all  entries  to  the  status  file  were  retained  for  30  days  at  the  inventory  control 
points  in  a concentrated  data  format.  A very  extensive  and  ambitious  programming 
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task  was  undertaken  at  FMSO  to  extract  these  data  and  generate  the  outputs  enumer- 
ated above.  A program  was  successfully  completed  and  a data  extraction  and  data 
analysis  effort  was  implemented  for  data  describing  30  days  of  operations  at  both 
SPCC  and  ESO.  The  data  analysis  was  segmented  with  respect  to  the  priority  of 
the  requisition  (Issue  Group  I,  n,  or  El),  the  extended  price  of  the  requisitioned 
items  (more  or  less  than  $2,200.00),  and  whether  it  was  an  FSN  versus  part  num- 
bered requisition.  A presentation  of  the  results  of  modeling  such  data  from  SPCC 
is  presented  in  Chapter  3. 

In  addition  to  satisfying  data  needs,  RMC  prepared  a detailed  analytic  model 
utilizing  a queueing  theoretic  approach  which  specified  elapsed  time  distributions 
for  requisitions  being  serviced  under  a multipriority  system.  Although  the  analytic 
principles  employed  were  general  to  any  number  of  priorities,  programs  were 
written  for  the  specific  cases  of  no  priorities,  two  priorities,  and  three  priorities. 

A discussion  of  how  priorities  were  taken  into  account  is  given  in  Chapter  3 with 
a detailed  description  of  the  mathematics  provided  in  Appendix  A. 

Other  portions  of  the  analysis  were  programmed  and  placed  in  on-line  status. 
Tests  were  provided  to  determine  if  the  work  unit  arrival  rate  distributions  and 
the  service  time  distributions  that  were  actually  observed  sufficiently  resembled 
those  distributions  which  were  assumed  the  case  in  order  to  parameterize  the 
multipriority  queueing  model.  Convolution  routines  were  programmed  to  assist 
in  combining  several  branch  throughput  time  distributions  into  a single  "path” 
throughput  time  distribution.  Finally,  automated  procedures  were  developed  which 
estimate  the  proportion  of  requisitions  leaving  a given  subdivision  that  eventually 
reach  each  other  subdivision  and  which  enumerate  the  various  "paths”  through  a 
facility  that  requisitions  might  take. 

CAPACITY  FACTORS 

An  input  to  the  multipriority  model  is  an  estimate  of  the  capacity  (work  units 
per  time  period)  of  each  functional  subdivision  of  a facility  to  process  requisitions. 
This  estimate  plus  the  assumed  work  load  are  used  by  the  multipriority  queueing 
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model  to  generate  throughput  time  distributions  for  each  such  subdivision.  In  the 
case  of  stock  points  the  MUACS*  data  base  provides  work  unit  standards  per  dedi- 
cated man-hour  for  many  tasks  performed  within  the  facility,  on  the  basis  of  engi- 
neering estimates.  Standards  for  tasks  for  which  engineering  data  do  not  exist  are 
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extrapolated  from  current  observations.  In  addition  the  DIMES  reporting  system 
provides  monthly  observations  of  work  units  completed  and  man-hours  expended 
which  can  readily  be  used  to  estimate  capacity  factors. 

The  data  situation  at  inventory  control  points  was  found  to  be  less  straight- 
forward. Since  the  resources  of  the  ICP  are  utilized  in  a number  of  functions 
other  than  requisition  processing,  performance  data  are  not  available  for  requisi- 
tion processing  ger  se.  However,  investigation  by  RMC  Research  Corporation 
personnel  revealed  that  the  cost  accounting  categories  maintained  at  SPCC  and 
ESO  would  permit  a fair  estimate  of  the  man-hours  dedicated  (charged)  to  requi- 
sition processing  within  each  functional  subdivision.  From  other  records  FMSO 
personnel  were  able  to  achieve  estimates  of  the  number  of  requisitions  processed 
by  each  subdivision.  In  combination,  these  data  were  used  to  estimate  man-hour/ 
capacity  relationships. 

OUTLINE  OF  THIS  REPORT 

Part  of  the  effort  to  implement  the  methodology  discussed  here  entailed  con- 
sulting support.  In  particular,  the  RMC  Research  Corporation  was  charged  with 
helping  FMSO  personnel  to  use  the  method  in  modeling  supply  system  operations 
other  than  ICP  requisition  processing  and  in  planning  and  designing  inputs  and  out- 
puts of  data  extraction  and  analysis  of  records  in  the  ICP  status  file.  The  perform- 
ance of  these  tasks  entailed  numerous  working  sessions  and  presentations  at  the 
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FMSO  facility  in  Mechanicsburg,  Pennsylvania  and  at  NAVSUP  headquarters. 


1.  Manpower  Utilization  and  Control  System. 

2.  Defense  Integrated  Management  Engineering  System. 

3.  Naval  Supply  Systems  Command. 
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This  report  is  organized  as  follows.  Chapter  2 presents  the  methodology 
utilized  and  some  Illustrations  of  its  application  to  data  describing  operations  at 
NSC  Norfolk.  Chapter  3 presents  analysis  of  a portion  of  the  data  collected  at 
SPCC  with  an  extrapolation  of  throughput  times  under  various  changes  in  capacity 
assumptions.  Chapter  4 presents  a discussion  of  the  insights  gained  from  the 
methodology  implementation.  Appendix  A is  a mathematical  discussion  of  the 
multipriority  queueing  model  and  other  statistical  procedures  developed  for  the 
methodology.  Appendix  B contains  listings  of  the  relevant  computer  programs. 
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METHODOLOGY 


The  method  chosen  to  study  the  resource  level/throughput  time  relationship  is  as 
follows : 

• deconsolidate  the  facility  under  study  into  subdivisions  consistent  with 
data  availability  and  the  various  tasks  performed  by  the  facility; 

• document  the  "paths"  a requisition  may  follow  among  these  subdivisions 
and  estimate  the  frequency  with  which  each  such  "path"  occurs; 

• for  each  subdivision  estimate  a distribution  of  throughput  times  taking 
into  account  available  resources  and  their  associated  capacity,  work 
load,  and  waiting  times  associated  with  administratively  determined 
policies; 

• amend  the  subdivision  throughput  time  distributions  as  dictated  by  sched- 
uling factors  associated  with  the  movement  of  requisitions  within  the 
facility  and  convolute  the  appropriate  subdivision  distributions  to  achieve 
a "path"  throughput  time  distribution;  and 

• take  the  average  of  the  "path"  throughput  time  distributions  weighted  by 
the  frequency  with  which  each  path  occurs  to  achieve  a facility  through- 
put time  distribution.  * 

Examining  the  steps  comprising  the  methodology  reveals  that  throughput  time 
is  sensitive  to 

• servicing  times, 

• waits  due  to  backlogs. 


1.  The  subject  matter  could  well  be  segmented  for  throughput  times  with 
respect  to  the  alternative  disposition  of  the  requisition  (e.g. , for  an  ICP:  the 
distribution  of  times  associated  with,  being  sent  to  a stock  point,  placed  on  back- 
order, canceled,  etc. ). 
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• pickup  and  delivery  schedules, 

• the  "paths"  the  requisitions  follow,  and 

• administratively  determined  waits. 


The  service  time  portion  of  throughput  time  is  not  at  issue  in  this  study.  That  is, 
the  service  rate  of  an  individual  is  not  under  study;  rather,  the  capacity  of  an 
administrative  subdivision  is  a function  of  the  number  of  man-hours  assumed  avail- 
able given  the  service  rate  per  man-hour.*  Waits  due  to  backlogs  are  (presumably) 
precisely  that  portion  of  throughput  time  which  is  to  be  explained  by  the  queueing 
theory  portion  of  the  method,  given  the  assumptions  about  capacities  and  work  load 
levels.  It  is  the  variation  in  waiting  times  due  to  backlogs  predicted  by  the  queue- 
ing theory  on  the  basis  of  variations  in  capacities  that  constitutes  the  central  re- 
source level/throughput  time  relationship  of  the  methodology.  In  practice,  the 
queueing  theory  predicted  very  modest  backlogs,  as  observed  capacities  were  very 
large  compared  to  work  loads.  Some  discussion  of  the  implications  of  these  results 
is  contained  in  Chapter  4.  Pickup  and  delivery  schedules  very  often  tend  to 
dominate  the  throughput  time  predicted  by  the  methodology.  If  capacities  are  large 
compared  to  work  loads,  then  the  queueing  theory  will  tend  to  predict  that  a requi- 
sition is  ready  to  go  on  the  next  scheduled  pickup  after  its  delivery  to  a servicing 
station.  Indeed  even  a cursory  view  of  supply  system  operations  leads  to  the  con- 
clusion that  requisition  processing  throughput  times  can  be  very  significantly  in- 
fluenced by  variations  in  messenger  service. 

The  various  "paths"  a requisition  may  follow  through  a facility  are  related  to 
resources  indirectly.  Obviously  if  requirements  to  service  a requisition  are 
amended  to  include  greater  or  fewer  servicing  units,  then  the  resources  utilized 
in  processing  a requisition  are  accordingly  greater  or  fewer.  If  the  number  of 
man-hours  dedicated  to  each  servicing  unit  remains  constant,  then  changing  the 
"paths"  utilized  or  the  frequency  with  which  a "path"  occurs  will  change  the  work 
load  requirement  for  the  servicing  units  involved.  Such  changes  will  impact  on 


1.  This  service  rate  might  be  from  an  engineering  standard  as  with  MUACS 
or  more  normally  imputed  from  observation. 
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throughput  time  by  implying  longer  or  shorter  waits  due  to  queues.  The  "paths" 
followed  by  requisitions  are  determined  by  policy  decisions.  As  a result,  experi- 
ments concerning  the  influence  on  throughput  time  of  changing  servicing  require- 
ments have  not  been  attempted.  However,  the  methodology  is  easily  able  to  support 
such  experiments  if  desired.  Administratively  determined  waits  (e.g. , waits  due 
to  the  competitive  requirements  of  the  procurement  process)  provide  additive 
constants  to  the  relevant  throughput  time  distributions.  Experiments  with  admin- 
istratively determined  waiting  times  are  beyond  the  scope  of  this  methodology 
i mple  mentation . 

STEPS  IN  THE  METHODOLOGY 

The  analysis  of  requisition  processing  accepts  the  following  data  inputs: 

• a definition  of  the  elements  or  subdivisions  of  the  system  under  study; 

• the  frequency  distribution  of  work  unit  flows  among  the  elements  of  the 
system,  e.g. , the  proportion  of  work  units  leaving  each  element  that 
immediately  go  to  each  other  element; 

• the  work  units  per  time  period  capacity  of  each  system's  element;  and 

• the  work  load  (in  each  priority  as  appropriate)  of  each  element. 

The  automated  analytical  portions  of  the  method  are  as  follows  in  the  order  of  their 
usage* 

Tests  for  Distributional  Characteristics 

The  multipriority  queueing  model  was  implemented  under  the  assumptions 
that  the  underlying  input  stream  for  each  priority  is  Poisson  and  that  service 
times  are  exponentially  distributed.  Procedures  have  been  automated  which 
evaluate  the  degree  to  which  observed  arrivals  and  service  times  are  described 
by  these  assumptions. 
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Analysis  of  Work  Load  Distributions 


A procedure  has  been  programmed  which  accepts  the  frequency  distribution  of 
flows  among  system  elements  and  extrapolates  them  to  predict  the  work  load  of 
each  systems  element  as  a percent  of  the  requisitions  leaving  any  other  given  sys- 
tems element.  The  gross  work  load  of  an  element  may  be  calculated  by  multiply- 
ing this  percent  times  the  number  of  requisition  arriving  in  the  system  at  each 
originating  element  and  adding  the  products  together. 

Estimation  of  Throughput  Time  Distributions  for  Each  System  Element 

The  multipriority  queueing  model  accepts  data  describing  average  work  loads 
for  each  priority  considered  and  capacity  for  each  element.  The  model  outputs  the 
expected  throughput  time  distribution  for  each  element  by  priority. 

Path  Analysis 

A procedure  has  been  automated  which  enumerates  the  possible  "paths"  a requi- 
sition may  follow  through  a facility  and  computes  the  frequency  with  which  each  path 
occurs.  This  analysis  is  based  on  the  data  describing  the  frequency  distribution 
of  requisition  flows  among  system's  elements. 

Convolution  Routines 

The  throughput  time  distribution  for  a "path"  followed  by  a requisition  is  deter- 
mined by  convoluting  the  throughput  time  distributions  for  each  of  the  constituent 
elements  of  the  path  with  adjustments  included  for  pickup  and  delivery  schedules  and 
administratively  determined  waiting  times. 

Using  the  data  and  automated  procedures  described  above,  the  steps  of  the  analy- 
sis described  below  are  as  follows: 

• eliminate  cycles  from  flow  matrix, 

• establish  work  load  factors, 

• establish  capacity  factors, 

• implement  multipriority  queueing  model, 

• path  analysis, 
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e establish  frequency  distribution  of  throughput  times  per  path,  and 
• calculate  distribution  of  system  throughput  times. 

STEP  I:  ELIMINATE  CYCLES  FROM  FLOW  MATRIX 

Data  input:  Frequency  distribution  of  flows  among  elements  (flow  matrix). 

Program:  Manual;  circumstances  such  that  a work  unit  flows  from  one  element 
to  another  and  returns  are  eliminated  by  treating  the  originating  element  dif- 
ferently (i.e.,  as  a distinct,  new  element)  upon  its  receipt  of  the  returning  work 
unit. 

Output:  Amended  flow  matrix  with  cycles  eliminated. 

Remark:  The  flows  among  elements  are  organized  in  a matrix  format.  If  a 
work  unit  is  sent  from  one  facility  element  to  another  (say  from  the  customer 
service  branch  to  the  technical  branch  of  a stock  point)  and  returns  a cycle  is 
established  which  the  automated  analytic  procedures  would  treat  endlessly. 

In  order  to  avoid  the  "error"  such  a process  would  introduce  into  the  analysis, 
"dummy"  elements  are  introduced  which  receive  requisitions  returning  to  an 
element  and  distribute  them  to  all  but  that  element  from  which  they  have  re- 
turned thus  eliminating  the  cycle. 

STEP  II:  ESTABLISH  WORK  LOAD  FACTORS 

Data  input:  Amended  frequency  distribution  of  flows  among  elements  (amended 
flow  matrix). 

Program:  Matrix  power  series. 

Output:  Work  units  reaching  any  given  element  as  a proportion  of  those  leaving 
any  given  element. 

Remark:  Of  primary  interest  is  the  work  load  per  element.  This  is  given  for 
each  element  by  taking  the  work  units  assumed  to  be  received  by  each  initial 
element  and  multiplying  by  the  proportion  of  those  receipts  which  were  received 
by  each  other  element.  However,  for  the  purposes  of  special  studies  the  pro- 
gram output  displays  receipts  by  an  element  with  respect  to  work  units  leaving 
any  other  element.  If  desired,  work  load  distinctions  can  be  maintained  with 
respect  to  a number  of  criteria.  As  discussed  in  Chapter  1,  the  data  describing 
SPCC  and  ESO  were  segmented  with  respect  to  issue  group,  price,  and  FSN 
versus  part  numbered  requisitions. 
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STEP  m:  ESTABLISH  CAPACITY  FACTORS 


The  determination  of  an  estimate  of  the  size  of  a system  subdivision  work 
unit  capacity  as  a function  of  dedicated  resources  is  difficult.  Several  approaches 
are  possible;  one  way  is  to  represent  rather  carefully  the  "technology"  or  process 
involved  at  a fine  level  of  detail  and  on  the  basis  of  the  identification  of  bottlenecks 
or  other  criteria  establish  capacities.  Such  a procedure  is  expensive  and  requires 
a level  of  detail  not  otherwise  usable  to  the  analysis.  The  DIMES  supplemental 
data  report  for  stock  point  operations  includes  "standards"  for  work  unit  comple- 
tions per  man-hours  expended.  Unfortunately  these  "standards"  have  not  been 
established  for  all  tasks  involved  in  stock  point  operations.  Further,  an  inspection 
of  the  data  reveals  that  actual  performance  is  often  very  different  (by  a factor  of 
100  percent  or  more)  than  that  which  would  be  predicted  by  an  application  of  the 
standards.  Individuals  on  the  spot  testified  that  such  variance  could  be  explained 
in  part  by  the  fact  that  many  of  the  standards  were  out  of  date  and  no  longer 
descriptive  of  current  practices.  In  any  event,  such  standards  were  not  available 
for  ICPs. 

An  alternative  procedure  is  to  observe  the  work  units  which  are  in  fact  proc- 
essed and  compare  this  figure  to  the  man-hours  dedicated  to  requisition  processing. 
Assuming  that  data  are  available  which  provide  this  information,  some  lower  bound 
to  capacity  can  be  inferred  (e.g. , what  is  observed  is  possible).  A complicating 
difficulty  in  the  assessment  of  capacity  is  that  of  "multi-use"  resources.  The  man- 
power resources  of  an  ICP  are  utilized  in  a number  of  activities  in  addition  to 
requisition  processing.  Cost  accounting  categories  reveal  (to  a degree)  the  man- 
hour split  among  various  activities  by  the  same  individual.  Of  principle  concern, 
however,  is  whether  or  not  a change  in  work  load  will  result  in  a change  in  waiting 
times  because  a given  man-hour  must  cope  with  more  or  fewer  requisitions,  or  will 
result  in  a change  in  man-hours  dedicated.  That  is,  will  more  requisitions  lead  to 
longer  waiting  times  due  to  queues  or  will  man-hours  be  released  from  other  tasks 
to  accommodate  the  increased  work  load  ? A number  of  studies  of  requisition 
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processing  have  recently  been  conducted.  A uniform  observation  of  these  studies 
is  that  peak  work  loads  are  much  larger  than  average  (or  "normal")  work  loads. 

The  general  conclusion  from  a study  of  these  observations  is  that  requisition  proc- 
essing resources  are  utilized  well  below  their  apparent  capacity  to  complete  work 
units.  When  such  low  levels  of  utilization  are  extended  into  models  which  Include 
an  analysis  of  queueing  discipline,  the  results  of  the  modeling  usually  fail  to  explain 
waits  due  to  backlogs  which  are  actually  observed.  Further  discussion  of  this  prob- 
lem will  be  offered  in  Chapter  4. 

For  purposes  of  illustrating  the  method,  data  describing  NSC  Norfolk  will  be 
presented.  FMSO  personnel  collected  cost  account  data  from  the  DIMES  supple- 
mental data  report  for  the  months  of  July,  August,  and  September  1972.  These 
data  permitted  the  calculation  of  work  units  per  man-hour  and  man-hours  per  calen- 
dar hour  (e.g.,  the  "apparent"  number  of  individuals  working  per  hour)  for  each 
system's  element  for  each  month.  The  measurement  of  capacity  is  calculated  by 
taking  the  product  of  the  maximum  observed  rate  of  processing  per  man-hour  and 
the  maximum  man-hour  allocations.  The  processing  rate  and  man-hour  allocations 
used  may  occur  in  different  time  periods.  Given  capacity  per  hour  and  given  the 

proportion  of  requisitions  received  by  NSC  Norfolk  that  reach  each  given  element 

2 

from  Table  2-1  (work  units  per  work  unit),“  the  maximum  number  of  requisitions 
received  by  NSC  Norfolk  per  hour  which  each  element  can  support  is  calculated  by 
dividing  hourly  capacity  by  the  work  units  per  work  unit  proportion  (e.g. , if  "pur- 
chase" can  pi*ocess  34.06  work  units  per  hour  an!  if  12  requisitions  per  thousand 
received  by  NSC  Norfolk  go  to  "purchase,"  then  "purchase"  can  support  a gross 
arrival  rate  of  34.06/.  012  = 2838.3  work  units  per  hour). 


1.  Some  data  sources  are:  the  SIMCOM  model  of  SPCC  developed  by  FMSO 

in  1971-72;  MUACS  data  collected  by  RMC  Research  Corporation  personnel  describ- 
ing operations  at  NSC  Norfolk  during  1971-72;  the  DIMES  supplemental  data  report 
for  NSC  Norfolk  for  July,  August,  and  September  1972;  SPCC  cost  accounting  data 
for  1972  and  SPCC  performance  data  for  January  and  February  1973;  and  Lynch  and 
Verich,  Requisition  Throughput  Time  Simulation  at  NSC  San  Diego,  March  1973. 

2.  e.g. , work  units  received  by  the  element  per  work  unit  received  by  the 
facility. 
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Table  2-1 


DATA  DESCRIBING  CAPACITY  FACTORS  FOR  NSC  NORFOLK 


Stock  Point 
Division 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(3) 

Processing 

Rate/ 

Hour 

Man-Hours/ 

Calendar 

Hours 

Capacity/ 

Hour 

Work  Units 
Reed  atDiv./ 
Work  Units 
Rec'd  at 
Stock  Point 

Maximum 
Supportable 
Work  Units/ 
Hour 

Relative 
Work  Load 

Observed 
Work  Load 
As  % of 
Capacity 

Work  Units 
Arriving 
Per 

Man-Hour 

Purchase 

.91 

37.43 

34.06 

.012 

• 2838.3 

.70 

.39 

.35 

Customer  Service 

22.22 

71.02 

1578.06 

.340 

4641.6 

.43 

.24 

5.33 

Technical 

2.81 

19.35 

54.37 

.020 

2718.5 

.73 

.41 

1.15 

Issue 

9.34 

189.57 

1770.58 

.705 

2511.5 

.79 

.44 

4.11 

Packing 

6.97 

177.16 

1234.80 

.620 

1991.6 

1.00 

.56 

3.90 

Going  through  this  calculation  for  each  element  reveals  substantial  differences 
in  the  gross  arrival  rates  that  each  element  can  support.  Of  the  five  elements 
modeled  by  the  queueing  formalism  the  maximum  supportable  arrival  (that  associ- 
ated with  ’’customer  service")  is  2.38  times  the  minimum  such  rate  (that  associated 
with  "packing").  As  a result,  "packing"  comprises  the  bottleneck  in  the  processing 
system  described  by  the  data.  Relative  to  its  capacity  "customer  service"  receives 
only  43  percent  of  what  "packing"  receives  per  hour.  These  relative  weightings  were 
used  as  appropriate  for  each  element. 

The  data  for  NSC  Norfolk  is  summaried  in  Table  2-1.  The  columns  of  Table 
2-1  are  as  follows: 

Column  1:  The  largest  of  the  three  numbers*  calculated  by  dividing  "work  units 
completed"  by  "man-hours  expended.  ” 

Column  2:  The  largest  of  the  three  numbers  calculated  by  dividing  "man-hours 
expended"  by  "man-hours  per  man  per  month  (i.e. , 8 x "work  days  in  month"). 

Column  3:  (Column  1)  x (Column  2). 

Column  4:  The  number  of  work  units  received  by  an  element  per  work  units  re- 
ceived by  NSC  Norfolk,  from  Column  1 of  Table  1. 

Column  5:  (Column  3)/ (Column  4). 

Column  6;  A rescaling  of  Column  5 such  that  each  entry  is  multiplied  by 
(1991. 6)_1.  Such  indicates  the  relationship  of  the  maximum  work  units  re- 
ceived to  that  allowed  by  the  bottleneck  capacity  of  "packing.” 

Column  7:  The  data  describing  operations  at  NSC  Norfolk  showed  "packing" 
receiving  an  average  of  56  percent  of  its  capacity.  The  implied  workload  of 
each  other  element  is  scaled  accordingly  (i.e.,  Column  7 = (Column  6)  x (.56)). 

Column  8:  The  steady  state  work  units  per  hour  per  hour  as  governed  by  "pack- 
ing" receiving  56  percent  of  capacity  = (Column  1)  x (Column  7). 

(N.B. , the  computer  and  the  "communications"  aspects  of  stock  point  operations 
were  not  modeled.  For  the  range  of  work  loads  to  be  considered,  It  was  determined 
that  requisitions  received  were  always  ready  to  move  on  per  the  pickup  and  delivery 
schedule. ) 


*i.e. , for  July,  August,  and  September. 
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STEP  IV:  IMPLEMENT  MULTIPRIORITY  QUEUEING  MODEL 


Data  input:  For  each  subdivision  of  the  system,  work  load  rates  per  time  period 
for  each  priority  from  STEP  n above;  element  capacities  from  STEP  m. 

Program:  Any  of  a no-priority,  two-priority,  and  three-priority  programs  are 
available.  Although  issue  group  distinctions  usually  provide  a context  of  three 
priorities,  it  is  sometimes  the  case  that  issue  groups  are  combined. 

Output:  For  each  priority, 

• the  average  number  of  work  units  in  the  system, 

• the  average  wait  in  the  system, 

o the  variance  of  the  system  wait,  and 
0 the  distribution  of  throughput  times. 

STEP  V : PATH  A NA  LYSIS 

Data  input:  Amended  frequency  distribution  of  flows  among  systems  elements. 

Program:  A path  analysis  program  traces  through  the  system  and  documents 
the  various  sequences  of  systems  subdivision  that  a requisition  can  encounter 
in  moving  through  the  system. 

Output:  For  each  terminating  event  (e.g. , manner  of  completing  the  requisi- 
tion), 

• an  enumeration  of  the  alternative  sequences  of  systems  subdivisions 
that  a requisition  can  encounter,  and 

0 the  frequency  with  which  each  "path"  occurs. 

Remark:  A path  is  a sequence  of  elements  that  a requisition  or  corresponding 
material  may  feasibly  encounter.  A path  is  always  initiated  by  an  element 
which  receives  the  work  unit  from  outside  the  system  (e.  g. , communications) 
and  terminated  by  an  element  which  passes  the  work  unit  on  outside  the  system 
or  otherwise  satisfies  the  requisition  (e.g.,  the  "customer  service"  subdivi- 
sion of  a stock  point  sends  a requisition  to  an  ICP  or  the  "purchase"  subdivi- 
sion of  an  ICP  places  a requisition  on  back  order).  Given  the  frequency  with 
which  a requisition  flows  between  all  pairs  of  systems  elements,  the  frequency 
with  which  a path  occurs  is  measured  by  taking  the  product  of  the  frequency 
values  of  the  flows  which  comprise  the  sequence. 
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STEP  VI:  ESTABLISH  FREQUENCY  DISTRIBUTION  OF  THROUGHPUT  TIMES  PER 
PATH 

Data  input:  Frequency  distribution  of  throughput  times  for  each  element  from 
Step  IV;  an  enumeration  of  path  elements  from  Step  V. 

Program:  Convolution  routine;  the  distribution  of  throughput  times  per  path  is 
computed  by  convoluting  the  throughput  time  distributions  of  the  path  elements; 
an  average  waiting  time  per  path  is  also  calculated. 

Output:  Frequency  distribution  of  throughput  times  and  average  throughput  time 
per  path. 

Remark:  Institutional  factors  with  respect  to  batching  or  other  system's  waits 
not  relating  to  the  actual  servicing  of  a work  unit  are  taken  into  account  at  this 
point.  The  waiting  times  per  work  unit  are  amended  to  account  for  such  waits 
as  appropriate . (e.g. , If  requisitions  leave  an  element  every  two  hours  and  the 
waiting  time  distribution  is  with  respect  to  each  hour,  then  an  amended  distribu- 
tion is  constructed  showing  zero  probability  on  the  odd  hours  and  the  sum  of  the 
current  and  the  hour  previous  probability  on  the  even  hours;  if  work  units  must 
wait  in  an  element  for  a fixed  average  time  independent  of  servicing,  then  the 
entire  waiting  time  distribution  is  "shifted"  forward  in  time  by  that  average 
wait. ) 

STEP  Vn:  CALCULATE  DISTRIBUTION  OF  SYSTEM  THROUGHPUT  TIMES 

Data  input:  Probability  of  path  occurrence  from  Step  V,  frequency  distribution 
of  throughput  times  per  path  from  Step  VI. 

Program:  Manual;  the  frequency  distributions  for  the  system  as  a whole  is  cal- 
culated by  taking  the  weighted  average  of  the  path  throughput  time  distributions 
with  the  path  probabilities  serving  as  weights.  In  addition,  the  frequency  distri- 
bution of  throughput  times  for  any  given  terminating  event  (e.g. , throughput 
times  in  a stock  point  for  requisitions  sent  to  an  ICP)  can  be  determined  by 
taking  the  weighted  average  of  the  relevant  paths  using  normalized  path  proba- 
bilities as  weights  (i.  e. , multiplying  the  relevant  path  probabilities  by  a factor 
such  that  they  sum  to  one). 

Output:  Frequency  distributions  of  throughput  times  and  average  throughput  time 
for  the  system  as  a whole  and  for  each  manner  in  which  a work  unit  leaves  the 
system  as  desired. 

Remark:  A waiting  time  frequency  distribution  can  of  course  be  calculated  in 
response  to  any  criterion  which  identifies  a subset  of  paths  in  the  system  (e.g. , 
all  paths  that  use  a given  element,  a given  two  elements,  etc. , . . . ). 
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ILLUSTRATION  OF  THE  METHODOLOGY 

Data  describing  requisition  processing  and  materials  handling  at  NSC  Norfolk 
were  collected  by  FMSO  personnel  and  will  be  used  here  to  illustrate  the  various 
techniques  and  procedures  of  the  methodology. 

STEP  I:  ELIMINATE  CYCLES 

Consideration  of  data  availability  and  the  nature  of  the  tasks  performed  at  a 
stock  point  suggested  that  stock  point  operations  should  be  represented  by  seven 
subdivisions  of  the  stock  point  "system.  " These  are 

• Communications, 

• Customer  Service, 

• "Computer, " 

• Technical, 

• Purchase, 

• Issue,  and 

• Packing. 

These  administrative  subdivisions  were  chosen  in  order  to  use  existing  cost  account 
data  as  reported  in  the  DIMES  supplemental  data  report  (i.e. , a finer  detail  is  not 
readily  accessible).  An  analysis  of  the  flows  among  these  elements  revealed  that 
requisitions  flowing  from  the  "computer"  to  "Technical"  would  in  part  flow  back  to 
the  "computer. " To  eliminate  this  cycle,  the  "computer"  was  represented  by  two 
elements:  "Computer  A"  which  accounts  for  requisitions  initially  received  and  "Com- 
puter B"  which  accounts  for  requisitions  returned  to  the  "computer"  from  "Technical. " 
Requisitions  returning  to  the  stock  point  from  DISC1  were  ignored  as  they  would  com- 
prise only  about  1 percent  of  the  total  (as  a result  total  throughput  time  for  the  stock 
point  as  a whole  will  fail  to  take  into  account  the  long  waiting  times  associated  with 


1.  Defense  Logistics  Services  Center. 
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filling  a requisition  cycling  through  DLSC;  however,  throughput  time  for  each  sepa- 
rate terminating  circumstance  is  unbiased).  The  following  flow  diagram,  Figure  2-1, 
summarizes  the  amended  flow  matrix. 

STEP  II:  ESTABLISH  WORK  LOAD  FACTORS 

Conceptually  there  are  a number  of  ways  to  determine  the  work  load  of  each  por- 
tion of  the  stock  point  as  a function  of  the  total  number  of  requisitions  arriving  at  the 
facility  per  same  time  period.  Simple  averaging  or  regression  procedures  suffer 
from  the  need  to  take  into  account  the  time  lags  between  a requisition  arrival  at  the 
stock  point  and  its  eventual  receipt  by  some  subdivision  of  the  stock  point.  The  pro- 
cedure used  here  is  that  of  tracing  through  the  network  illustrated  in  Figure  2-1  and 
determining  for  each  system's  element  the  proportion  of  work  units  it  receives. 

This  would  be  a very  tedious  process  even  for  only  moderately  large  systems. 

As  a result  a computer  program  was  prepared  which  computes  the  proportion 
of  requisitions  reaching  any  given  systems  element  relative  to  those  leaving  each 
other  systems  element.  This  program  utilizes  the  frequency  distributions  of  work 
unit  flows  between  pairs  of  systems  elements  such  as  those  displayed  in  Figure  2-1. 
For  a stock  point  all  requisitions  are  represented  as  arriving  at  the  "communications” 
subdivision.  As  a result  the  work  load  of  each  other  subdivision  may  be  calculated 
by  determining  the  proportion  of  requisitions  leaving  "communications"  which  reaches 
each  other  subdivision  and  multiplying  that  proportion  by  the  number  of  requisitions 
assumed  to  be  arriving  at  "communications."  The  calculation  of  proportions  is  dis- 
played in  Table  2-2  for  NSC  Norfolk  based  upon  the  frequency  distributions  of  flows 
given  in  Figure  2-1.  Generally,  the  numbers  in  the  table  indicate  the  number  of 
requisitions  reaching  the  row  coordinate  as  a proportion  of  the  number  leaving  the 
column  coordinate.  As  "communications"  is  always  encountered  first,  the  first  column 
of  Table  2-2  contains  the  proportions  of  particular  interest.  For  example,  the  first 
entry  reveals  that  34  percent  of  the  requisitions  received  by  "communications"  are 
sent  on  to  "customer  service,"  and  so  on.  It  should  be  noted  that  the  analysis  takes 
the  entire  network  into  account  and  includes  all  possible  paths  between  "communica- 
tions" and  each  other  subdivision  when  computing  the  relevant  proportion. 
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Communications 


Figure  2-1:  NETWORK  OF  REQUISITION  AND  MATERIAL  FLOWS  AT  NSC  NORFOLK 
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Table  2-2 


FLOW  PROPORTIONS  AMONG  STOCK  POINT  DIVISIONS 


CLSSRV  CCVU  _ 

C. 3400 

CCVPA CCVU  CLSSRV  

1.0000  I. 0000 

J £C  H Q C V 0 C LSSRV  CC  VP  A 

C . 0200  0.02C0  C.C200 

CCVPB  _ CCVL_.__COS.SkV CO  VP  A JF:CH 

0.CC60  0.0060  0.0060  0.3000 

I S SUE_  _CCVU CL S SR V _CC_VPA TECH COVPB _ 

C.705L  6 . 7 05 1 0.705 l 6.2550  0.0500 

PACK  _ CCfU CUSSRV CCVPA  TFCH COMPP  1SSUF  

C .6205  0.6205  0.6205  0.2244  0.7480  0.8800 

CISC  __  _CCV,U  __ CUSSRV CCVPA 

C.OICO  0.0100  C.OIOO 

PLRCH  CCVU  _ CL SSRV  CCVPA  _TECH 

0.0124  0.0124  0.0124  0.6200 

ICP  _ CCVU CL  SSRV CP  MPA  TECH CGVPe 

C .2725  o’.27?5  C.2725  0.1250  0.1500 

SHIP  CCVU  CLSSKV  CCVPA  TFCH  _ COVPP.  ISSUC  PACK 

C.7C51  C.7051  C.7051  0.2550’  0.8500  1.0000  1.0000 


STEP  m:  ESTABLISH  CAPACITIES 

As  discussed  above  an  accurate  estimate  of  an  element's  capacity  based  on 
available  secondary  data  is  difficult  to  achieve.  Even  when  engineering  estimates 
are  available  they  are  often  incomplete  or  out  of  date.  As  a result,  it  is  necessary 
to  estimate  capacity  factors  from  observations  of  performance.  Presumably  the 
capacity  per  man-hour  of  an  element  is  no  smaller  than  observed  work  rates.  In 
Table  2-1,  work  unit  capacities  per  man-hour  are  provided  based  on  the  maximum 
per  hour  work  rate  and  maximum  quantity  of  man-hours  per  calendar  hour  expended 
over  a three-month  period.  In  the  multipriority  queueing  model  output  displayed 
below  (starting  on  page  2-17)  for  the  issue  and  packing  elements,  the  capacities  are 
given  respectively  by: 

capacity  for  issue  = 9.34  per  man-hour,  and 
capacity  for  packing  = 6.97  per  man-hour.  * 


1.  These  measurements  were  derived  from  data  collected  for  July,  August, 
and  September  1972. 
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STEP  IV:  IMPLEMENT  MU LTIPRIORITY  QUEUEING  MODEL1 

An  issue  not  firmly  resolved  by  this  study  is  the  choice  of  a unit  of  analysis  for 
the  queueing  theoretic  model.  That  is,  what  is  to  comprise  the  basic  servicing 
unit  (e.g. , how  many  man-hours  per  calendar  hour  comprise  a "servicing  unit  hour")? 
A study  of  intra-element  structure  was  deemed  beyond  the  scope  of  the  analysis  and  no 
other  statistical  means  of  approaching  this  issue  was  developed.  Generally,  the 
choice  of  a unit  of  analysis  was  governed  by  the  goal  of  fitting  the  output  of  the  multi- 
priority  queueing  model  to  observed  behavior.  Further  discussion  of  this  problem  is 
provided  in  connection  with  the  modeling  of  SPCC  operations  discussed  in  Chapter  3. 

The  multipriority  queueing  model  was  applied  to  the  NSC  Norfolk  data  for  the 
elements  "issue"  and  "packing"  using  the  man-hour  as  the  basic  unit  of  analysis.  This 
choice  is  recommended  by  the  fact  that  issue  and  packing  operations  have  a very 
simple  structure  similar  to  many  parallel  servicing  units.  The  capacity  and  arrival 
rate  data  were  not  further  scaled  as  comprehensive  performance  data  were  not  avail- 
able for  the  stock  point  at  the  time  the  example  was  formulated.  Using  the  capacity 
estimates  given  above  and  the  arrival  rates  given  above  in  Table  2-1  the  queueing 
model  was  implemented  for  three  cases: 

© the  base  case  per  the  observed  data,  and 

• two  "experimental  cases  assuming  a 20  and  a 40  percent  increase  in 
work  load. 

The  arrivals  were  split  among  the  priorities  as  follows: 

• IG-I,  10  percent; 

• IG-n,  40  percent;  and 

• IG-m,  50  percent. 

The  outputs  of  the  model  are  self-explanatory.  Statistics  for  the  several  priorities, 
high  to  low,  are  displayed  from  left  to  right,  or  from  top  to  bottom  as  appropriate. 

1.  The  multipriority  queueing  model  is  given  in  Appendix  A. 
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STEPS  V AND  VI:  PATH  ANALYSIS  AND  PATH  FREQUENCY  DISTRIBUTIONS 


For  purposes  of  example  only  the  single  sequence  of  events  leading  to  the  issue 
and  packing  of  stocked  material  in  response  to  stock  numbered  requisitions  was  con- 
sidered. Referring  to  the  diagram  given  in  Figure  2-1,  this  sequence  entails  the  re- 
ceipt of  the  requisition  by  the  communications  element;  sending  the  requisition  to  the 
computer  whether  or  not  via  the  customer  service  element;  and  then  sending  the 
requisition  to  the  issue  element  and  the  corresponding  material  to  the  packing  ele- 
ment. Hie  general  procedure  is  to  estimate  the  throughput  time  distribution  for  the 
sequence  of  events  by  convoluting  the  throughput  time  distributions  for  each  element 
in  the  sequence. 

For  the  issue  sequence  in  this  illustration  the  elements  "communications,” 
"customer  service, " and  "computer  A"  were  consolidated  into  a single  element. 
Generally  it  was  supposed  that  a requisition  would  be  processed  by  each  of  these 
elements  and  sent  on  in  a fashion  governed  by  pickup  and  delivery  schedules.  For 
the  sake  of  the  example  these  were  assumed  to  be 

IG-I:  once  per  hour, 

IG-n:  twice  a day,  and 
IG-ni:  once  a day. 1 

The  distributions  were  then  convoluted  with  the  throughput  time  distributions  for 
issue  and  packing  for  each  Issue  (group  and  for  each  of  the  three  cases:  base  case, 
+20  percent  arrivals,  and  +40  percent  arrivals.  These  throughput  time  distributions 
are  displayed  on  the  following  pages . 


1.  In  fact,  batching  requisitions  for  the  warehouse  occurs  once  a day  for  both 
IG-n  and  IG-ni  although  at  different  times.  IG-I  requisitions  are  not  sent  quite  as  often  as 
once  every  work  day  hour.  In  addition,  IG-I  requisitions  are  batched  several  times 
during  the  night. 
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Inspection  of  this  illustration  reveals  that  even  an  increase  of  40  percent  in  work 
load  leads  to  a projection  by  the  model  of  a very  small  increase  in  throughput  time. 

As  mentioned  above,  this  conclusion  is  typical  of  studies  of  this  kind  and  can  be  attri- 
buted to  the  (apparently)  low  levels  of  utilization  at  supply  system  activities.  In  the 
analysis  of  SPCC  data  provided  in  the  next  chapter  more  comprehensive  performance 
data  were  available  and  the  data  were  scaled  such  that  the  multipriority  queueing 
model  more  closely  fitted  observed  throughput  time  distributions.  In  these  instances 
the  projections  by  the  model  were  more  sensitive  to  changes  in  the  capacity-arrival 
rate  relationship. 

STEP  VII:  SYSTEM  THROUGHPUT  TIMES 

If  more  than  one  path  is  studied,  the  throughput  time  distribution  for  the  system 
is  achieved  by  taking  the  weighted  average  of  the  path  distributions.  The  weight 
associated  with  a path  is  the  frequency  with  which  the  path  occurs  as  calculated  by 
the  product  of  the  frequency  of  the  individual  flows  which  make  up  the  path.  For  the 
issue  sequence  studied  the  path  frequency  would  be  .617  (e.g. , 61.7  percent  of  the 
requisitions  received  by  NSC  Norfolk  follow  this  path)  as  calculated  by  the  flow 
frequencies  given  in  Figure  2-1. 
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MODELING  REQUISITION  PROCESSING  AT  SPCC 

The  steps  in  the  modeling  procedure  presented  in  Chapter  2 were  arranged  in  the 
sequence  in  which  they  arose  in  the  model  building  process.  This  chapter  reports 
on  the  implementation  of  the  model.  As  it  turned  out  the  model  was  applied  with 
the  steps  occurring  in  a somewhat  different  sequence.  As  a result  the  steps  are 
presented  in  this  chapter  in  the  sequence  of  their  application.  The  numbering  of 
the  steps  presented  in  Chapter  2 is  retained. 

The  application  of  the  methodology  to  inventory  management  practices  was  handi- 
capped by  a lack  of  supporting  secondary  data.  The  general  requirement  was  a docu- 
mentation of  where  requisitions  flowed  within  an  organization,  how  long  they  remained 
at  each  of  the  various  elements  which  processed  them,  and  an  accounting  of  the  man- 
hours expended  for  requisition  processing.  No  currently  collected  data  base  explicitly 
responded  to  these  needs.  RMC  Research  Corporation  personnel  with  close  support 
and  cooperation  from  personnel  of  the  Administrative  Management  Division  of  SPCC 
determined  that  the  currently  maintained  "requisition  status  file"  at  SPCC  appeared 
to  contain  data  sufficient  to  support  the  analysis,  though  its  reporting  format  would 
require  some  amendment  to  generate  precisely  the  data  required. 

Generally,  it  was  desired  to  use  the  status  file  to  document  the  "histories"  of 
requisitions  within  SPCC  both  in  terms  of  where  the  requisitions  flowed  and  how  long 
they  remained  at  the  various  subdivisions  of  SPCC  which  processed  them.  It  was 
also  hoped  that  man-hour  allocations  to  requisition  processing  could  be  determined 
at  the  same  time  from  cost  accounting  data  maintained  separately.  The  difficulty  in 


using  the  status  file  stemmed  from  the  fact  that  each  time  the  file  was  updated  with 
respect  to  a requisition's  location  or  status,  all  information  about  its  previous  history 
was  erased.  Some  400  requisitions  were  selected  and  their  descriptions  in  the  status 
file  were  extracted  once  per  day  for  a 30-day  period.  As  a result  of  this  exercise  it 
was  determined  that  the  data  contained  in  the  status  file  were  otherwise  sufficient  to 
support  the  analysis  if  some  means  could  be  found  to  preserve  them. 

Accordingly,  FMSO  personnel  studied  the  problem  of  using  the  status  file  at 
SPCC  to  generate  the  needed  requisition  histories.  It  was  determined  that  all 
entries  to  the  status  file  were  maintained  in  a concentrated  format  data  file  for 
roughly  four  weeks.  A very  substantial  computer  programming  task  was  under- 
taken by  FMSO  personnel  to  extract  the  needed  data  from  the  concentrated  file  and 
perform  the  required  data  analysis  to  generate  the  inputs  for  the  methodology.  The 
extraction  and  analysis  programs  were  successfully  applied  to  files  maintained  at 
both  SPCC  and  ESO  for  data  describing  a 30-day  period  documenting  all  requisi- 
tions being  processed  by  the  inventory  control  points.  The  data  analysis  program 
genei'ated  the  following  descriptors  of  ICP  operations: 

• for  each  subdivision  of  the  facility  the  number  of  requisitions  received 
per  day  over  the  period; 

• for  each  subdivision  the  number  of  requisitions  completed  on  each  day; 

• the  proportion  (and  number)  of  requisitions  which  flowed  directly  from 
each  subdivision  to  each  other  subdivision; 

• for  each  subdivision  the  frequency  distributions  of  final  actions,  e.g. , 
the  proportionate  breakdown  of  what  eventually  happened  to  requisitions 
which  reached  each  given  subdivision; 

• the  distributions  of  total  elapsed  times  required  to  complete  requisitions 
broken  down  with  respect  to  type  of  final  action;  and 

• the  distributions  of  times  requisitions  spent  in  each  subdivision. 

The  data  were  segmented  with  respect  to  each  of  the  three  issue  groups  and  with 
respect  to  unit  prices  above  and  below  $2, 200. 00.  Only  requisitions  not  referred 
automatically  were  included. 
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STEP  I:  ELIMINATE  CYCLES  FROM  FLOW  MATRIX 


The  local  routing  codes  used  in  the  SPCC  status  file  enabled  the  study  of  opera- 
tions at  SPCC  at  the  branch  level.  Data  extraction  and  analysis  were  performed  for 
the  following  branches : 


Code 

(Division)  Branch 

494 

"Computer" 

(720) 

(Financial  Control) 

722 

Material  Accounting 

(770) 

(Purchase) 

771 

Selected  Items  Purchasing 

773 

Buying 

774 

Contract  Management 

778 

Purchase  Services 

(810) 

(Support  Determination) 

813 

Ordnance  Support 

814 

Technical  Support 

(840) 

(Stock  Control) 

841 

Special  Support 

842 

Customer  Requirements 

844 

Safety  and  Electrical 

845 

Machinery  and  Engine 

846 

Weapons 

(870) 

(Nuclear  Equipment  Support) 

871 

Nuclear  Weapons 

872 

Nuclear  Propulsion 

873 

Nuclear  Propulsion  Material 

(890) 

(Strategic  Sys terns  Support) 

891 

Support  Determination 

892 

Material  Management 

893 

Program  Management 

To  reduce  the  scope  of  the  data  modeling  effort  it  was  decided  to  limit  the  modeling 
effort  to  federal  stock  numbered  requisitions  with  unit  cost  below  $2, 200.  An  exami- 
nation of  requisition  flows  among  branches  revealed  that  a requisition  could  flow 
through  a given  branch  a second  time  in  a number  of  instances.  As  a result  a number 
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of  "dummy"  branches  were  introduced  as  discussed  above.  Figure  3-1  displays  the 
pattern  of  flows  among  branches  for  federal  stock  numbered  requisitions  with  unit 
price  below  $2,200  at  SPCC  as  revealed  by  the  data  extraction  and  analysis  program. 
The  "dummy"  branches  are  designated  by  adding  a letter  to  the  branch  code.  Branch 
81  is  a consolidation  of  branches  813  and  814,  84  a consolidation  of  844,  845,  and  846, 
and  89  a consolidation  of  891,  892,  and  893. 

The  flow  matrix,  Figure  3-1,  was  then  analyzed  to  document  the  total  flows. 
Computer  software  was  developed  which  determined  the  proportion  of  requisitions 
leaving  any  given  branch  which  arrive  at  each  other  branch.  This  proportion  and 
the  number  of  requisitions  entering  the  system  can  be  combined  to  provide  a 
"steady  state"  estimate  of  work  load  per  branch.  However,  since  work  load  data 
were  collected  from  the  status  file,  this  steady  state  estimate  was  not  used  in  the 
multipriority  queueing  model.  Table  3-1  summarizes  the  results  of  this  flow 
analysis  for  the  various  final  actions: 

back  order  (BKORD), 
sent  to  ICP  (ICP), 
referred  (REFRD), 
purchased  (SPOTBUY), 
cancelled  (CANC),  and 
a substitute  found  (SUB). 

The  final  actions  are  the  row  coordinates  and  the  branch  codes  the  column  coordi- 
nates. Each  entry  gives  the  proportion  of  requisitions  received  by  a branch  that  is 
completed  by  each  of  the  actions  indicated.  Since  essentially  all  requisitions  are 
processed  through  the  computer  the  entries  for  branch  "494"  display  the  distribution 
of  final  actions  for  the  ICP  as  a whole.  Below  this  display  is  an  enumeration  of  the 
proportion  of  requisitions  received  by  SPCC  (and  not  referred  automatically)  that 
goes  to  each  branch. 

STEP  V:  PATH  ANALYSIS 

Since  the  number  of  branches  which  could  potentially  be  studied  is  large,  a 
path  analysis  was  conducted  next  to  identify  those  branches  appearing  prominently 
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Figure  3-1:  SPCC:  FLOW  OF  REQUISITIONS,  IG-H,  FSN,  UNIT  PRICE  < $2, 200 


I 


Table  3-1 

PROPORTION  OF  REQUISITIONS  PROCESSED  PER  BRANCH  PER  FINAL  ACTION 


BKOMO 

49* 

722 

771 

773 

773A 

774 

77*A 

778 

778A 

81 

8 1 A 

8* 

84  A 

0.127* 

0.1700 

0.0657 

0.1128 

0.1*05 

0.1838 

0.2500 

0.0865 

0.109* 

0.1522 

0.2275 

0.5350 

0.6500 

Bk0»0 

8*2 

04?  A 

B72 

872A 

673 

873A 

89 

* 

0.2*6* 

0.1300 

0.2215 

0.3600 

0.1616 

0.2000 

0.3534 

icp 

*94 

771 

773 

77  3A 

77* 

77*A 

778 

778A 

01 

81 A 

8* 

842 

842A 

0.0196 

0.0226 

0.0133 

0.0229 

0.0199 

0.0500 

0.0166 

0.0272 

0.0986 

0.1400 

0.0087 

0.1096 

0.1600 

1CP 

872 

873 

89 

0.0126 

0.0054 

0.1197 

REFRO 

49* 

722 

771 

771 A 

773 

773A 

77* 

774  A 

778 

778A 

81 

81 A 

8* 

0.1618 

0.0612 

0.1546 

0.1*00 

0.1681 

0.1667 

0.1770 

0.1500 

0.1521 

0.1658 

0.1585 

0.1939 

0.2331 

MEFRO 

84  A 

6*2 

8*2A 

872 

872a 

873 

873A 

89 

0. 3500 

0.2990 

0.5600 

0.32*1 

O.S500 

0.1567 

0.2700 

0.0*07 

SPOTPUY 

494 

722 

771 

77 1 A 

773 

773A 

77* 

774A 

778 

778A 

61 

81 A 

8* 

0.6346 

0.2*86 

0.7059 

0.6600 

0.6*19 

0.603* 

0.6162 

0.5500 

0.7121 

0.693* 

0.5828 

0.4386 

0.2219 

SPOTPUY 

0*2 

8*2  A 

872 

673 

873A 

89 

0.3127 

0.1000 

0.3929 

0.4623 

0.5300 

0.1053 

CANC 

49* 

722 

771 

773 

773A 

770 

8*2 

6*2A 

672 

873 

89 

• 

0.0037 

0.0700 

0.00*9 

0.0069 

0.0056 

0.0026 

0.0035 

0.0500 

0.0063 

0.0287 

0.00*0 

SUB 

773A 

77* 

778A 

8* 

672 

672A 

0.0020 

0.0031 

0.00*2 

0.0013 

0*0018 

0.0700 

PROPORTION  OF  TOTAL  REQUISITIONS  PROCESSED  BY  EACH  BRANCH 


494 

722 

771 

77 1 A 

773 

773A 

77* 

774A 

778 

77dA 

81  « 

611  64  84* 

0.0672 

0.5200 

0.29*6 

0.2000 

0.0953 

0.0*20 

0.0292 

0.3000 

0.167* 

0.2020 

0.0337  0.1624  0.0259 

494 

8*2 

072 

872A 

873 

073A 

89 

BK0R0 

ICP 

REFRO 

SPOTBUY 

CANC 

0.0312 

0.0150 

0.0066 

0.0183 

0.0066 

0.0022 

0.1274 

0.0196 

* 0.1618 

0.63*8 

0.0037 

- 


I*  «k.-  - 


in  requisition  processing  histories.  The  methodology  includes  a computer  routine 
which  utilizes  the  flow  matrix  to  enumerate  the  various  possible  paths  and  computes 
the  frequency  with  which  each  path  occurs.  Inspection  of  Figure  3-1  reveals  that 
the  number  of  all  possible  paths,  however  infrequent,  would  be  very  large.  As  a 
result  the  routine  was  amended  to  enumerate  only  those  paths  which  occur  at  least 
once  per  thousand  requisitions.  These  paths  were  organized  by  final  action  and  are 
displayed  at  the  end  of  the  chapter.  Paths  occurring  at  least  once  per  thousand  requi- 
sitions accounted  for  83  percent  of  all  requisitions.  If  only  paths  occuring  once  per 
hundred  are  counted,  55  percent  of  all  requisitions  are  accounted  for.  The  very 
large  number  of  paths  generated  by  the  flow  matrix  indicates  that  the  level  of  detail 
achieved  using  the  branches  of  SPCC  as  a unit  of  analysis  is  somewhat  finer  than 
was  originally  supposed.*  Paths  which  occur  at  least  once  per  hundred  requisitions 
are  enumerated  in  Table  3-2.  It  was  decided  to  go  forward  only  with  an  analysis  of 
branches  appearing  in  the  paths  given  in  Table  3-2.  These  are  771,  773,  778,  81, 
and  84. 

STEP  m:  ESTABLISH  CAPACITY  FACTORS 

Since  work  load  estimates  were  taken  directly  from  the  data  extraction  and 
analysis  efforts,  rather  than  inferred  from  the  flow  analysis,  observed  man-hour 
allocations  were  used  to  determine  work  load  per  (potentially)  available  man-hour. 

As  a result  it  is  appropriate  to  consider  performance  data  on  which  to  base  capacity 
estimates.  Cost  accounting  data  for  the  year  1972  were  examined  to  document 
branch  performance.  For  each  month  two  figures  were  calculated  (when  possible): 

• the  average  number  of  requisitions  completed  per  man-hour  charged  to 
requisition  processing,  and 

• the  average  number  of  man-hours  charged  to  requisition  processing  per 
calendar  hour. 


1.  There  are  104  paths  occurring  once  per  thousand  requisitions  which  fail 
to  account  for  17  percent  of  the  requisitions  processed  (of  those  not  automatically 
referred  by  the  computer).  A more  complete  application  of  the  method  to  so  large 
a system  would  require  automation  of  additional  aspects  of  the  path  analysis:  a 
straightforward  task. 
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Table  3-2 


REQUISITION  PROCESSING  PATHS 
OCCURRING  AT  LEAST  ONCE  PER  HUNDRED  REQUISITIONS 


mmmmamammmm 

Path  No. 

Frequency 

Branches 

Final  Action  E 

1 

.050 

494,  84 

Back  Order  I 

i 2 

.042 

494,  778,  771A 

Referred 

3 

.012 

494,  773 

Referred 

4 

.018 

494,  84 

Referred 

5 

.229 

494,  771 

Spot  Buy 

6 

.031 

494,  778,  771A 

Spot  Buy  | 

7 

.032 

494,  771,  778,  771A 

Spot  Buy 

8 

.030 

494,  773,  778,  771A 

Spot  Buy 

9 

.013 

494,  84,  778A,  771A 

Spot  Buy 

' 10 

.029 

494,  81,  771A 

Spot  Buy  | 

11 

.028 

494,  771,  81,  771A 

Spot  Buy 

12 

.044 

494,  773 

Spot  Buy 
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As  it  happened  these  calculations  were  performed  for  all  branches  and  the  results  are 
presented  in  Table  3-3  below.  The  servicing  rate  per  maivhour  of  a branch  is 
chosen  as  the  maximum  observed  rate  in  Table  3-3  and  the  available  man-hours  (e.g. , 
man-hour  capacity)  per  calendar  hour  are  chosen  as  the  maximum  observed  alloca- 
tion in  Table  3-3. 

STEP  H:  ESTABLISH  WORK  LOAD  FACTORS 

Arrivals  at  each  branch  in  each  priority  were  provided  by  the  data  extraction 
and  analysis  program.  It  was  provisionally  decided  to  express  arrivals  in  terms  of 
the  (potentially)  available  man-hours  as  determined  from  the  19?2  performance  data 
given  in  Table  3-3.  This  was  achieved  by  dividing  monthly  arrivals  by  an  assumed 
176  calendar  hours  in  the  month  and  dividing  the  quotient  by  the  maximum 
observed  man-hour  allocation  for  1972.  As  discussed,  this  procedure  treats  the 
branch  intra-structure  as  essentially  independent.  The  results  of  these  calculations 
are  provided  in  Table  3-4. 

STEP  IV:  IMPLEMENT  MULTIPRIORITY  QUEUEING  MODEL 

At  this  stage  in  the  analysis  the  multipriority  queueing  model  is  applied.  The 
model  accepts  parameter  values  measuring  the  average  arrival  rate  and  servicing 
rate  per  time  unit  for  each  of  three  priorities  and  outputs  throughput  time  distribu- 
tions and  other  statistics  as  displayed  in  Chapter  2.  In  order  to  implement  the  model 
some  assumption  must  be  made  about  branch  capacities  (as  discussed)  and  the  "time 
unit  of  analysis"  must  be  selected  (e.g. , arrivals  and  capacity  Per  hour?  per  day? 

. . .).  The  goal  was  to  implement  the  model  using  parameters  descriptive  of  the 
time  period  viewed  by  the  data  extraction  and  analysis  effort  and  compare  the  model 
results  with  observed  throughput  time  distributions.  The  observed  throughput  time 
distributions  revealed  average  times  of  several  days  magnitudo  and  throughput  time 
distributions  covering  over  10,  sometimes  over  20,  days. 

As  discussed  it  was  provisionally  decided  to  model  the  branch  man-hour,  that 
is,  express  capacities  and  arrivals  in  terms  of  the  projected  available  man-hours 
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Table  3-3 


PROCESSING  RATES  AND  MAN-HOUR  ALLOCATIONS  PER  BRANCH  FOR  1972 


S Branches 

Jan. 

Feb. 

Mar. 

Apr. 

May 

June 

July 

Aug. 

Sept. 

Oct. 

Nov. 

Dec.  H 

S 771 

1.93 

36.99 

36.12 

3.06 

41.82 

1.00 

39.64 

2.56 

39.60 

1.68 

37.61 

1.88 

33.44 

2.01 

34.20 

1.64 

37.53 

1.71 

34.82 

1.80 

33.29 

1.76  B 
29.38  1 

773 

.37 

77.10 

.35 

78.84 

.36 

89.67 

.31 

84.29 

.32 

73.06 

.31 

75.20 

.36 

73.64 

.43 

79.61 

.40 

78.66 

.38 

76.29 

.19 

82.00 

1 

.20  | 
75.69  B 

774 

.31 

58.02 

.48 

56.57 

.54 

66.49 

.45 

58.75 

.49 

52.90 

.52 

52.57 

.51 

49.52 

.64 

53.22 

.56 

54.10 

.46 

52.77 

.36 

52.47 

.33  | 

44.26  I 

778 

.63 

64.94 

.73 

62.18 

.60 

71.72 

.46 

73.39 

.68 

66.22 

.62 

60.06 

.70 

59.16 

.81 

60.00 

.67 

60.77 

.72 

59.35 

.87 

60.73 

.69  | 

52.21  I 

813 

1.10 

20.79 

.73 

19.14 

1.23 

16.65 

1.16 

16.81 

1.24 

14.01 

1.32 

15.48 

1.62 

16.89 

1.25 

15.26 

1.55 

17.53 

1.51 

14.61 

1.67 

14.29 

1. 78  I 
11.87  3 

814 

.78 

75.23 

.81 

70.05 

.80 

68.97 

1.12 

61.22 

.71 

55.95 

.74 

50.81 

.71 

54.75 

.84 

55.19 

.72 

55.95 

.73 

55.67 

.70 

58.43 

.74 

47.47 

841 

.92 

2.79 

.96 

2.05 

. 66 
1.89 

.98 

4.94 

3.99 

5.73 

4.40 

4.08 

5.41 

2.42 

6.64 

3.64 

5.53 

3.02 

5.52 

2.58 

6.84 

3. 18  | 
6.14 

842 

2.64 

57.05 

3.53 

54.65 

3.79 

61.63 

2.96 

56.60 

4.18 

52.65 

4.94 

49.41 

2.83 

47.21 

2.39 

54.64 

4.37 

53.65 

2.40 

47.92 

aTMan 

4.63 

51.56 

4.54 

48.19 

Format: 


requisitions  per  man-hour 


man-hours  per  calendar  hour 
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Table  3-3  - Processing  Rates  and  Man-Hour  Allocations  Per  Branch  for  1972  (Continued) 


Branches 

Jan. 

Feb. 

Mar. 

Apr. 

May 

June 

July 

Aug. 

Sept. 

Oct. 

Nov. 

1 i 

Dec.  a 

844 

2.38 

2.  70 

2.57 

3.76 

4.05 

3.82 

3.93 

4.25 

3.43 

5.67 

2.85 

4.32 

26.82 

24.37 

25.68 

23.55 

21.47 

18.06 

15.98 

16.33 

16.96 

17.52 

20.08 

19.90  j 

845 

1.39 

1.66 

1.47 

3.48 

2.87 

1.83 

1.94 

1.86 

1.40 

2.85 

1.51 

1.45 

23.96 

22.73 

26.06 

22.67 

20.77 

19.20 

20.20 

20.74 

23.43 

21.31 

21.21 

18.30  : 

846 

2.28 

3.42 

— 

3.79 

4.20 

3.93 

3.60 

3.03 

— 

3.15 

3.91 

3.26 

31.01 

29.25 

30.02 

24.92 

22.86 

20.72 

22.08 

21.27 

24.20 

23.38 

23.83 

21.29 

871 

1.07 

1.  60 

1.71 

1.58 

1.65 

— 

1. 26 

1.24 

1.10 

1.24 

1.44 

1.87 

1.58 

- 

1.95 

1.50 

0.95 

0.58 

0.65 

0.45 

0.62 

0.65 

0.44 

872 

1.41 

1.64 

— 

1.90 

1.72 

2.01 

1.97 

1.76 

1.45 

2.39 

1.85 

1.61 

10.58 

10.54 

- 

9.98 

10.61 

7.83 

8.39 

8.96 

9.24 

7.55 

7.90 

7.38  ; 

873 

— 

— 

.. 

— 

0.13 

0.24 

0.15 

0.33  j 

5.38 

5.28 

- 

5.41 

5.30 

5.26 

5.00 

5.01 

5.84 

6.35 

7.01 

6.40  I 

891 

3.03 

3.35 

3.06 

2.54 

3.21 

4.73 

2.37 

2.41 

3.58 

3.45 

2.09  J 

12.34 

12.86 

11.15 

9.77 

8.41 

9.52 

11.09 

9.57 

8.80 

9.23 

8.94 

8.26  « 

892 

1.26 

20.23 

0.93 

24.74 

1.04 

23.35 

1.06 

22.07 

0.95 

19.09 

1.11 

19.05 

1.02 

22.52 

0.92 

20.32 

1.06 

19.69 

1.23 

19.73 

0.76  I 
16.38  1 

893 

«_ 

_ 

- 

- 

- 

- 

- 

0.40 

0.05 

- 

- 

- 

“ 1 

Format: 


requisitions  per  man-hour 


man-hours  per  calendar  hour 
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Table  3-4 


OBSERVED  REQUISITION  ARRIVALS  PER  MONTH  AND  REQUISITION  ARRIVALS  PER  AVAILABLE  MAN-HOUR 


Branch 

771 

773 

778 

Issue  Group 

I 

n 

in 

1 

n 

m 

I 

n 

m 

i 

n 

ni 

Arrivals  per  month 

2,318 

4,218 

1,697 

671 

1,779 

4,573 

1,607 

3,803 

3,145 

262 

586 

773 

Arrivals  per  branch 
man-hour  available 

.315 

.573 

.231 

.043 

.113 

.290 

.124 

.294 

.243 

.072 

.160 

.211 

Branch 

814 

844 

845 

846 

Issue  Group 

I 

n 

m 

I 

n 

m 

I 

n 

HI 

I 

n 

m 

Arrivals  per  month 

942 

1,386 

1,553 

318 

1,350 

841 

208 

767 

394 

258 

1,199 

974 

Arrivals  per  branch 
man-hour  available 

.071 

.105 

.117 

.067 

.286 

.178 

.045 

.167 

.086 

.047 

.220 

.178 

per  calendar  hour.  The  best  observed  servicing  rate  in  1972  was  chosen  as  the 
capacity  measure  and  the  maximum  (monthly)  allocation  during  that  period  as  the 
measure  of  available  man-hours.  Servicing  rates  were  set  equal  for  all  priorities. 
Arrivals  per  man-hour  were  set  at  those  average  rates  observed  (in  January, 
February  1973)  using  the  data  extraction  and  analysis  program.  The  results  were 
disappointing  though  not  entirely  unexpected.  Expressed  in  hourly  terms  the 
multipriority  queueing  model  depicted  throughput  times  of  only  a few  hours.  In 
no  instance  did  the  model  generate  distributions  of  times  which  stretched  over 
many  days  as  did  the  observed  data.  These  results  confirmed  earlier  findings, 
even  though  now  the  queueing  model  was  additionally  taking  into  account  the  priority 
system  created  by  the  three  issue  groups. 

When  the  inability  to  model  long  throughput  times  was  first  encountered  in  the 
methodology  development  phase  of  the  study,  deficiencies  in  model  assumptions  were 
immediately  suspected.  Accordingly,  statistical  procedures  were  found  which  tested 
the  assumptions  of  a Poisson  input  stream  and  exponentially  distributed  servicing 
times  (it  was  not  in  any  case  supposed  that  discrepancies  in  distributional  assump- 
tions alone  could  explain  the  very  large  waiting  times  observed).  These  procedures 
are  discussed  in  Appendix  A.  As  it  turned  out  servicing  times  (as  opposed  to 
elapsed  times  in  a branch)  were  not  directly  observed  and  the  test  of  their  distri- 
butional characteristics  was  not  conducted.  However,  arrivals  were  documented 
and  the  sizes  of  daily  arrivals  were  tested  against  daily  arrival  sizes  projected 
under  the  assumption  that  they  were  Poisson  distributed.  The  Kolmogorov- Smirnov 
test  was  employed.  The  distributions  of  arrivals  were  tested  at  all  but  two  branches 
at  SPCC  (data  problems  eliminated  two  branches  which  were  not  in  any  case  modeled). 
At  the  5-percent  level  of  significance  the  Poisson  assumption  was  not  rejected  for 
any  of  the  branches  tested. 

Next  the  capacity  estimates  employed  were  called  into  question.  Generally,  the 
branches  under  study  conduct  a number  of  activities  in  addition  to  requisition  proc- 
essing. Perhaps  the  use  of  the  best  observed  performance  entailed,  practically 
speaking,  an  over-estimate  of  servicing  ability.  In  the  spirit  of  using  the  mathe- 
matical formalism  descriptively,  it  was  decided  to  search  for  capacity  assumptions 
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which  would  give  descriptive  results  when  incorporated  in  the  mathematical  model. 
Further,  it  was  realized  that  the  special  handling  procedures  for  IG-I  requisitions 
implied  a smaller  capacity  than  that  for  IG-H  and  IG-m  requisitions.  The  refine- 
ment of  alternative  capacities  was  also  incorporated  into  the  model.  The  effort  to 
fit  capacities  failed  to  produce  realistic  results  and  was  abandoned.  In  general  the 
problem  was  the  need  to  specify  extremely  high  processing  rates  for  IG-n  and 
IG-in  requisitions  in  order  to  compensate  for  the  lower  capacity  for  IG-I  requisi- 
tions . 

Although  it  was  known  that  many  of  the  waiting  times  experienced  in  requisition 
processing  are  not  directly  related  to  available  processing  resources,  it  was  diffi- 
cult to  accept  that  the  independence  between  resources  and  elapsed  times  was  so 
extreme  that  the  multipriority  queueing  model  was  totally  unable  to  document  any, 
even  tenuous,  connection.  As  a result,  a final  fitting  of  the  model  was  implemented. 
If  the  servicing  and  arrival  rates  are  scaled  down  vuth  respect  to  shorter  time  units, 
the  throughput  time  distributions  depicted  by  the  model  become  fuller  and  extend  over 
more  (smaller)  periods.  This  characteristic  of  the  multipriority  queueing  model  was 
used  to  replicate  observed  throughput  time  distributions  as  follows.  The  IG-II  and 
IG-m  servicing  rate  was  set  at  the  maximum  observed  rate  for  1972.  The  IG-I 
servicing  rate  was  set  at  one-half  that  value.  The  observed  arrivals  were  expressed 
per  available  branch  man-hour  per  calendar  hour  using  the  maximum  man-hour 
allocation  observed  in  1972.  The  output  of  the  multipriority  queueing  model  was 
relabeled  in  terms  of  days.  The  observed  throughput  time  distribution  for  IG-n 
requisitions  was  studied  and  the  number  days  required  for  90  percent  of  the  requi- 
sitions to  be  completed  was  noted.  A scalar  was  then  sought  such  that  if  the  arrival 
and  servicing  rates  of  the  queueing  model  were  multiplicatively  scaled  by  that  scalar, 
the  IG-n  throughput  time  distribution  would  predict  all  requisitions  completed  in  the 
number  of  days  within  which  90  percent  was  actually  observed  to  be  completed.  For 
purposes  of  comparison  the  observed  distribution  was  then  normalized  (e.g. , re- 
scaled proportionately  to  show  100  percent  rather  than  90  percent  completed  in  the 
observed  number  of  days). 
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Several  remarks  are  immediately  appropriate.  The  degree  to  which  a mathe- 
matical model  generated  descriptively  can  then  be  used  predictively  is  an  empirical 
issue  for  which  this  study  was  not  designed  and  which  it  could  not  afford  to  investi- 
gate if  such  an  investigation  was  to  be  expensive  in  time.  It  was  hoped  that  the 
mathematical  model  could  depend  in  large  part  on  structural  (causal)  insights  and  a 
sufficient  verification  be  achieved  from  data  collected  over  several  months.  No 
use  of  the  model  could  be  so  verified.  Further,  the  procedure  to  "fit"  the  model 
to  observations  outlined  above  was  chosen  because  it  was  inexpensive  in  time  and 
other  resources.  The  effort  to  fit  the  model  was  undertaken  at  a late  hour  in  the 
contract  schedule.  It  would  have  been  no  less  (and  no  doubt  more)  appropriate  to 
scale  the  data  towards  making  the  observed  and  modeled  average  throughput  times 
as  nearly  equal  as  possible.1  Alternatively  the  absolute  deviation  between  the 
observed  and  modeled  distributions  could  provide  a measure  of  "fit. " The 
Kilmogorov-Smirnov  test  statistic  was  computed  for  several  of  the  branch  dis- 
tributions. An  alternative  fitting  procedure  would  be  the  reach  for  that  scale 
factor  which  minimized  that  statistic. 

Stated  simply  the  analytic  thrust  of  the  study  was  to  model  requisition  proc- 
essing based  on  a fair  estimate  of  the  characteristics  of  that  process  with  the 
resulting  advantage  of  an  a priori  focus  on  the  appropriate  mathematical  struc- 
tures. Observed  data  failed  to  verify  the  usefulness  of  the  model  so  constructed. 
The  further  contention  that  a substantial  portion  of  requisition  processing  is  in 
fact  modeled,  but  that  the  model  (for  whatever  reason)  requires  a certain  adjust- 
ing to  accurately  describe  reality,  goes  beyond  what  the  study  resources  had  been 
allocated  to  deal  with.  Hence,  the  following  use  of  the  model  is  not  offered  as  the 
best  means  of  fitting  the  model  to  the  available  data. 


1.  It  would  have  taken  more  time  however. 
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OUTPUTS  OF  THE  MULTIPRIORITY  QUEUEING  MODEL 

Using  the  data  and  procedures  outlined  above  the  multipriority  queueing  model 
was  run.  The  outputs  are  presented  in  the  following  tables.  Displays  are  given 
for  branches  771,  773,  778,  and  81  (the  weighted  average  of  branches  813  and 
814).  In  addition  to  the  base  case  described  above,  available  capacities  were 
changed  by  +10  percent  and  -10  percent  and  the  model  was  rerun.  There  were 
two  exceptions.  Branch  778  capacities  could  not  be  scaled  down  by  as  much  as 
10  percent  so  the  run  for  reduced  capacity  was  made  at  95  percent  base  case 
capacity.  In  addition,  Branch  773  achieved  a work  load  performance  which  ex- 
ceeded the  maximums  observed  in  1972.  In  this  instance  the  assumed  capacities 
were  scaled  up  to  those  (relatively)  of  the  most  highly  utilized  branch  (778).  As 
a result,  the  reduced  capacity  run  for  this  branch  was  also  made  at  95  percent. 
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Branch  771,  IG-I 


Branch  773,  IG-I 


Days 

Probabilities 

Observed 

Model: 
Base  Case 

Model: 

+10  Percent  Capacity 

Model  | 

-10  Percent  Capacity^ 

| 1 

.219 

.105 

.105 

.084 

2 

.212 

.105 

.105 

.084 

3 

.049 

.105 

.105 

.084  | 

4 

.024 

.105 

.167 

.084  I 

5 

.042 

.061 

.061 

.144 

6 

.103 

.061 

.061 

.059  | 

7 

.079 

.061 

.061 

.059 

8 

.054 

.061 

.128 

.059 

8 9 

.042 

.128 

.085 

.059 

10 

.012 

.069 

.036 

.103 

11 

.006 

.052 

.030 

.058 

12 

.012 

.024 

.019 

.036 

13 

.073 

.017 

.010 

. 024 

14 

.018 

.012 

.010 

.017 

15 

.037 

.009 

.005 

.012 

16 

.012 

.007 

.006 

.009 

17 

.006 

.005 

.003 

.007 

18 

.006 

.003 

.005 

19 

.003 

.004 

20 

.003 

.003 

.21 

.003 

22 

23 

24 

B 25 

26 

B Average  Days 

5.46 

6.28 

5.64 

6.66 
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Branch  778,  IG-I 


I 
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Branch  81,  IG-I 
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Branch  771,  IG-n 


Probabilities 

Days 

Observed 

Model: 
Base  Case 

Model: 

+10  Percent  Capacity 

Model  6 

-10  Percent  Capacity! 

1 

.159 

.141 

.211 

.141 

2 

.091 

' .141 

.211 

.141 

3 

.034 

.265 

.186 

.141 

4 

.053 

.124 

.186 

.099  j 

5 

.048 

.124 

.111 

.099 

6 

.104 

.085 

.046 

.099 

8 7 

.198 

.053 

.020 

.102 

8 

.142 

.028 

.014 

.072 

9 

.046 

.014 

.009 

.031 

10 

.020 

.012 

.006 

.026 

11 

.023 

.006 

.013 

12 

.039 

.006 

.013 

13 

.018 

.003 

.007 

14 

.025 

.007 

15 

.004 

16 

.003 

17 

.003 

18 

19 

20 

.21 

22 

23 

24 

25 

26 



||  Average  Days 

5.96 

3.85 

3.16 

12 J 
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Branch  773,  IG-II 


Days 

Probabilities  8 

Observed 

Model: 
Base  Case 

Model: 

+10  Percent  Capacity 

Model  B 

-10  Percent  CapacityB 

8 1 

.226 

.105 

.105 

.084 

' 2 

.169 

.105 

.105 

.084 

i 3 

.041 

.105 

.105 

.084 

I 4 

.050 

.105 

.180 

.084 

5 

.052 

.044 

.074 

.120 

6 

.054 

.044 

.074 

.035 

7 

.127 

.044 

.074 

.035 

8 

.034 

.044 

.103 

.035 

9 

.028 

.161 

.058 

.035 

10 

.019 

.084 

.036 

.161 

11 

.045 

.035 

.030 

.061 

12 

.015 

.036 

.015 

.058 

13 

.021 

.024 

.011 

.026 

14 

.039 

.017 

.010 

.027 

15 

.021 

.008 

.006 

.013 

16 

.004 

.010 

.005 

.015 

17 

.025 

.007 

.004 

.008 

18 

.009 

.006 

.004 

.009 

19 

.009 

. 003 

.005 

20 

.011 

.004 

.006 

.21 

.003 

.003 

22 

.002 

.003 

23 

.003 

24 

.002 

25 

.002 

26  .... 

■ Average  Days 

5.93 

6.64 

5.49 

i 
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Branch  778,  IG-n 


Days 

Probabilities  1 

Observed 

Model: 
Base  Case 

Model: 

+10  Percent  Capacity 

Model  fl 

-10  Percent  Capacity  h 

1 

.345 

.141 

.211 

.141 

2 

.268 

.141 

.211 

.141 

3 

.024 

.222 

.148 

.141 

4 

.026 

.081 

.148 

.081 

5 

.040 

.081 

.161 

.081 

6 

.027 

.155 

.052 

.081 

7 

.033 

.072 

.032 

.128  | 

8 

.035 

.039 

.015 

.085 

9 

.025 

.023 

.011 

.030 

10 

.030 

.015 

.007 

.030  | 

11 

.020 

.013 

.003 

.015  I 

12 

.031 

.007 

.014 

13 

.038 

.005 

.007  | 

14 

.028 

.004 

.007 

15 

.030 

.001 

.004 

16 

.004 

17 

.003 

18 

19 

20 

.21 

22 

23 

24 

25 

26 

. 

• 

1 Average  Days 

4.38 

4.29 

3.34 

— 
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Branch  81,  IG-II 
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The  Kolmogorov-Smimov  test  statistic  was  computed  for  each  base  case  against 
the  observed  distribution.  The  results  are  given  below. 

Table  3-5 

KOLMOGOROV-SMIRNOV  TEST  STATISTIC 


The  indicated  fits  are  not  too  good;  however,  there  are  several  extenuating  consider- 
ations. Hie  mathematical  model  included  a smoothing  routine  which  would  show 
badly  in  a test  of  this  kind  against  more  irregular  distributions.  Some  of  these  irregu- 
larities would  be  neutralized  in  the  convoluted  distributions  so  that  it  may  be  expected 
that  the  path  throughput  time  distributions  "fit"  the  model  better. 

STEP  VI:  ESTABLISH  FREQUENCY  DISTRIBUTIONS  OF  THROUGHPUT  TIMES  PER 
PATH 

The  convolution  of  the  branch  throughput  time  distributions  for  a number  of  paths 
are  displayed  below.  The  weighted  averages  of  path  distributions  were  not  taken. 
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A 


* 


Path:  494,  778,  771A,  Referred 
Frequency:  . 0416 


1 

Probabilities 

Days 

1G 

-I 

IG-n 

Observed 

Model: 

Observed 

Model: 

Model: 

Model: 

j 

Base 

Base 

+10  Percent 

-10  Percent 

‘ Average 

■ Days 

6.28 

5.70 

10.34 

8.12 

6.59 

9.77 

* 

►.  1 ’ 

2 

.114 

.044 

.055 

.020 

.045 

.020 

• 

l 3 

.197 

.125 

.074 

.040 

.089 

.040 

ij  4 

.111 

.155 

.040 

.089 

.115 

.060 

V 

V 

5 

.071 

.178 

.034 

.097 

.141 

.065 

* 

{ 6 

.069 

.171 

.040 

.116 

.154 

.071 

* 

7 

.073 

.130 

.059 

.111 

.133 

.076 

u 

• 

t 8 

.095 

.086 

.108 

.110 

.108 

.091 

9 

.085 

.053 

.117 

.107 

.082 

.096 

10 

.049 

.029 

.073 

.081 

.052 

.088 

\ 11 

.026 

.016 

.043 

.065 

.035 

.076 

$ 

12 

.025 

.008 

.039 

.049 

.023 

.064 

13 

.031 

.003 

.046 

.036 

.013 

.055 

14 

.027 

.001 

.046 

.026 

.008 

.048 

j 

f 

15 

.011 

.043 

.018 

.004 

.039 

> 

t 16 

.003 

.035 

.013 

.002 

.029 

1 17 

.003 

.024 

.008 

.022 

f 

18 

.002 

.021 

.005 

.016 

19 

.001 

.021 

.003 

.012 

i 

20 

.021 

.002 

.009 

21 

.019 

.001 

.006 

22 

.015 

.005 

*! 

23 

.009 

.003 

24 

.005 

.002 

| 25 

.003 

.001 

« 

! 26 

- .003 

• 

{ 

27 

.002 

- 

$ 

28 

.001 

t 

29 

*' 

f 

30 

?] 

f 31 

; 

5 

V 

\ 32 

* 

!> 

i L3 

% 

j 

34 

35 

- 

i 36 

| 37 

f 

t 

38 

t 

•j 

39 

) 

40 

ei 
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Path:  494,  773,  778,  771A,  Spotbuy 
Frequency:  .0297 


1 

Probabilities 

Days 

IG 

-I 

1 

IG-II 

Observed 

Model: 

Observed 

Model: 

Model: 

Model: 

* 

Base 

Base 

+10  Percent 

-10  Percent 

. Average 

• Days 

11.74 

11.98 

16.26 

14.68 

12.12 

17.10 

2 

0 

0 

0 

0 

0 

0 

3 

.026 

.005 

.013 

.002 

.005 

.002 

4 

.068 

.018 

.023 

.006 

.019 

.005 

t 

5 

.072 

.034 

.024 

.016 

.026 

.010 

6 

a- 

.051 

.053 

.020 

.026 

.044 

.015 

•. 

t 

t* 

7 

.045 

.069 

.029 

.037 

.062 

.022 

8 

.057 

.077 

.030 

.046 

.076 

.028 

t 

9 

.075 

.079 

.051 

.052 

.085 

.035 

10 

.080 

.077 

.064 

.057 

.092 

.040 

e 

H 

.070 

.075 

.057 

.061 

.092 

.044 

> 12 

.054 

.077 

.044 

.064 

.088 

.050 

13 

.044 

.076 

.045 

.069 

.082 

.057 

* 

i w 

.044 

.074 

.042 

.069 

.073 

.057 

■* 

ii  15 

.052 

.067 

.056 

.068 

.062 

.059 

16 

.051 

.057 

.054 

.065 

.051 

.059 

> 

U 17 

.035 

.045 

.048 

.061 

.041 

.057 

18 

.033 

.033 

.048 

.056 

.031 

.057 

'l 

19 

.027 

.024 

.039 

.048 

.024 

.055 

| 

20 

.023 

.017 

.039 

.040 

.017 

.051 

: 

21 

.021 

.013 

.037 

.033 

.013 

.046 

22 

.017 

.009 

.036 

.027 

.009 

.040 

I 

23 

.017 

.007 

.033 

.022 

.007 

.031 

« 

i 

r 24 

.010 

.005 

.029 

.017 

.005 

.031 

i 

25 

.007 

.003 

.022 

.013 

.003 

.027 

i 

t 

f 26 

.005 

.002 

' .021 

.010 

.002 

.022 

! 

27 

.004 

.001 

.019 

.008 

.001 

.018 

* 

* 

! 28 

.004 

.017 

.006 

.015 

> 

; 29 

.002 

.015 

.004 

.012 

V 

30 

.001 

.011 

.003 

.009 

31 

.009 

.002 

.007 

t 

32 

.007 

.002 

.006 

i 

33 

.006 

.001 

.005 

i. 

f 

34 

.005 

.003 

l 

l 

35 

.004 

.003 

p 

36 

.003 

.002 

4 

37 

.002 

.001 

l 

M 

38 

.002 

i 

t 

39 

.001 

r 

40 

.001 

f; 

•< 
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Path:  494,  771,  778,  771A,  Spotbuy 
Frequency:  .0310 


lh'obabilities 

• Days 

l 

IG 

-r 

IG-II 

Observed 

Model: 

Observed 

Model: 

Model: 

Model: 

t 

t 

Base 

Base 

+10  Percent 

-10  Percent 

Average 

r Days 

10.52 

8.83 

16.30 

11.99 

9.64 

14.58 

S 2 

0 

0 

0 

0 

0 

0 

i 

i 3 

.030 

.009 

.009 

.003 

.009 

.003 

! 4 

.077 

.036 

.017 

.008 

.028 

.008 

K ^ 

.074 

.067 

.015 

.023 

.051 

.017 

» 

G 

.063 

.102 

.014 

.039 

.079 

.022 

* 

! 7 

.059 

.131 

.017 

.061 

.105 

.033 

CO 

.064 

.142 

.025 

.075 

.120 

.042 

!,  9 

.081 

.135 

.042 

.089 

.116 

.052 

10 

.099 

.117 

.061 

.097 

.100 

.062 

l 11 

.075 

.091 

.055 

.098 

.080 

.069 

* 

: 12 

.058 

.065 

.047 

.095 

.061 

.070 

: 13 

.058 

.043 

.047 

.085 

.044 

.073 

■ 

\ 14 

.053 

.027 

.057 

.074 

.030 

.077 

1 15 

.054 

.016 

.068 

.061 

.020 

.070 

! 16 

.047 

.009 

.071 

.050 

.013 

.065 

a 

l 17 

.030 

.005 

.062 

.039 

.008 

.059 

■ 

18 

.021 

.002 

.048 

.030 

.005 

.052 

T 

\ 19 

.018 

.001 

.042 

.022 

.003 

.045 

* 

20 

.015 

.042 

.016 

.003 

.037 

j 21 

.013 

.041 

.012 

.001 

.031 

1 

22 

.008 

.038 

.008 

.026 

* 

! 23 

.004 

.032 

.005 

.020 

*' 

24 

.003 

.026 

.004 

.016 

! 

25 

.002 

.021 

.002 

.013 

\ 26 

.002 

' .019 

.001 

.010 

* 

l 27 

.001 

.017 

.007 

£ 

\ 28 

.015 

.005 

i 

! 29 

.012 

.004 

[ 30 

.009 

.003 

. 

f 31 

.006 

.002 

■ 

5 

\ 32 

.005 

.001 

V 

. 

| 33 

.004 

i 

i 

34 

t 35 

.003 

.002 

i 

t 

* 

i 36 

.001 

* 

37 

r 

2^L 

CO 

00 

9 

; 

39 

; 

i 40 

1— J|m 

f 
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Path:  494,  771,  81,  771A,  Spotbuy 
Frequency:  . 0282 


Probabilities 

Days 

r 

i 

•> 

IG 

-I 

IG-II 

Observed 

Model: 

Base 

Observed 

Model: 

Base 

Model: 
+10  Percent 

Model: 
-10  Percent 

J 

Average 
( Days 

11.44 

10.90 

17.11 

11.33 

9.40 

13.66 

i 2 

0 

0 

0 

0 

0 

0 

i 3 

.015 

.006 

.003 

.004 

.008 

.003 

‘ 

4 

.040 

.019 

.008 

.011 

.025 

.008 

2 5 

.060 

.036 

.010 

.027 

.050 

.017 

*; 

6 

.070 

.056 

.011 

.044 

.078 

.026 

: 

7 

.062 

.076 

.014 

.063 

.102 

.034 

* 

i 8 

.062 

.090 

.020 

.077 

.117 

.042 

t 

9 

.074 

.099 

.031 

.091 

.113 

.056 

1 

10 

.080 

.102 

.042 

.097 

.117 

.063 

1 11 

.073 

.099 

.046 

.098 

.095 

.070 

i I2 

.071 

.092 

.046 

.091 

.072 

.073 

13 

.065 

.080 

.051 

.082 

.058 

.072 

2 

| 14 

.054 

.066 

.059 

.070 

.047 

.071 

\ 

; 15 

.052 

.052 

.066 

.057 

.030 

.068 

1 

' 16 

.041 

.039 

.070 

.045 

.020 

.062 

l 

1 

’ 17 

.036 

.028 

.067 

.034 

.014 

.055 

1 

' 18 

.030 

.019 

.061 

.025 

.009 

.048 

! 19 

.023 

.013 

.058 

.019 

.005 

.041 

1 20 

.019 

.009 

.056 

.013 

.003 

.034 

j)  21 

.017 

.006 

.052 

.009 

.002 

.028 

22 

.012 

.003 

.046 

.006 

.001 

.022 

23 

.008 

.002 

.038 

.004 

.017 

Is 

24 

.005 

.001 

.031 

.003 

.013 

V 

25 

.004 

.026 

.002 

.010 

. 

26 

.003 

' .022 

.008 

7 

‘ 27 

.002 

.018 

.006 

| 28 

.001 

.014 

.004 

j 29 

.010 

.003 

t 

j 30 

.007 

.002 

1 31 
32 

.005 

.003 

.002 

.001 

i 

r 

j 33 

! 34 

.003 

.002 

i 

\ 

\ 35 

.001 

1 

36 

37 

% 

38 

r 

39 

40 

t 

l 
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Path:  494,  81,  771A,  Spotbuy 
Frequency:  . 0288 


Days 

I 

Probabilities 

IG 

-I 

IG-II 

Observed 

Model: 

Base 

Observed 

Model: 

Base 

Model: 
+10  Percent 

Model: 
-10  Percent 

Average 

j 

Days 

7.21 

7.78 

11.15 

7.56 

6.34 

9.04 

2 

.069 

.030 

.021 

.028 

.040 

.020 

i 

3 

.097 

.060 

.040 

.054 

.080 

.040 

4 

.134 

.086 

.032 

.091 

.120 

.062 

5 

.131 

.103 

.035 

.098 

.142 

.066 

* 

i 

6 

.077 

.111 

.041 

.108 

.145 

.071 

) 

7 

.087 

.109 

.051 

.114 

.128 

.082 

8 

.089 

.110 

.074 

.118 

.100 

.094 

\ 

t; 

9 

.084 

.098 

.094 

.097 

.073 

.094 

10 

.070 

.083 

.083 

.076 

.050 

.084 

11 

.049 

.067 

.069 

.059 

.034 

.071 

T 

12 

.034 

.049 

.070 

.044 

.021 

.060 

13 

.032 

.032 

.071 

.030 

.014 

.058 

J 

14 

.028 

.021 

.067 

.021 

.009 

.043 

* 

15 

.022 

.015 

.056 

.014 

.005 

.032 

J 

16 

.011 

.010 

.042 

.009 

.003 

.024 

17 

.006 

.006 

.032 

.006 

.001 

.018 

18 

.004 

.004 

.028 

.004 

.013 

*! 

19 

.003 

.002 

.025 

.002 

.010 

> 

20 

.001 

.001 

.019 

.001 

.007 

3 

21 

.012 

.004 

L 

22 

.009 

.003 

r 

% 

i 

23 

.005 

.002 

> 

: 

J 

24 

.004 

.002 

25 

.002 

f 

i 

26 

' .002 

i 

i 

27 

F 

V 

28 

t 

t. 

29 

J 

t 

30 

r; 

\ 

31 

1 

32 

1 

33 

l 

34 

i 

35 

l 

36 

l 

37 

t 

i 

38 

* 

»* 

y 

39 

i 

40 

i 

3-30 


PATH  ENUMERATION 

The  following  is  the  enumeration  of  paths  occurring  more  than  once  per  thousand 
requisitions.  The  paths  are  displayed  as  follows: 

• the  final  action  is  given  as  "user  ending  node"; 

• all  paths  start  at  "494,  " the  computer; 

• the  path  is  then  enumerated  backwards  starting  with  the  "user  ending  node,  " 
the  branches  are  shown  pairwise  followed  by  the  frequency  of  the  flow  be- 
tween the  pair;  and 

• the  frequency  of  the  path  is  displayed  as  "freq.  prod."  which  is  the  product 
of  the  pairwise  frequencies. 
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Table  3-G 


AN  ENUMERATION  OF  REQUISITION  PROCESSING  PATER 


USSR  ENDING 

■ NOOE 

IS: ♦ 22 

9K0RD 

starting  node 

IS: 

1 494 

NO. 

OF  pairs: 

3 

2d  d<0R0 

2 722 

0.170  2 

722 

3 

771 

0.070 

3 771 

1 494 

0.520 

freo.  ppod.  = 

0.0061880 

» , 

starting  node 

IS:  * 

1 494 

NO. 

OF  PAIRS: 

3 \ * V * 

22  6K0R0 

2 722 

0.170  2 

72? 

5 

773 

0.060 

5 773 

1 494 

0.200 

. 

FREO.  PROO.  = 

0.0027200 

STARTING  NODE 

is: 

1 494 

*N0. 

OF  PAIRS: 

4 

22  BKORO 

8 774A 

0.250  8 

774A 

10 

778A 

0.180 

10  776 A 

13  84 

0.320  13 

84 

1 

494 

0.100 

FREO.  PROD,  = 

0. 0014<.00 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS: 

2 

22  BKORD 

13  84 

0.500  13 

84 

1 

494 

0.100 

FREO.  PROD,  - 

0.0500000 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

3 

22  e^ORO 

13  84 

0.500  13 

84 

5 

773 

0.090 

5 773 

1 494 

0.200 

FREC.  PROD.  * 

« 

0.0090000 

STARTING  NODE 

is: 

1 4Q4 

NO. 

OF  PAIRS: 

3 

22  BKORD 

13  84 

0 .500  13 

84 

11 

61 

0.220 

11  81 

1 494 

0.080 

FR£C.  PROD.  * 

0.0088000 

STAwTJsg  node 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

22  B*OPD 

13  84 

0.500  13 

94  - * 

11 

61 

0.220 

11  81 

3 771 

0.150  3 

771 

1 

494 

0.520 

FREO.  PROD.  = 

0.0085800 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS*. 

4 

22  B<OPD 

. 13  84 

0.500  13 

84 

11 

81 

0.2  20 

1)  81 

5 773 

0.070  5 

773 

1 

494 

f 0.200 

FREO.  PRCD..S 

0.0015400 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

22  B<OPD 

13  *84 

0.500  13 

84 

11 

81 

0.220 

1 1 81 

9 778 

0.100  9 

770 

1 

494 

0.100 

FREO.  PROD.  * 

0.0011000 

starting  nooe 

IS: 

1 494 

NO. 

OF  PAIRS: 

5 

22  B«ORD 

13  84 

0.500  13 

84 

11 

61 

0.220 

1 t ei 

9 77d 

0.100  9 

778 

3 

771 

0.200 

3 771 

1 494 

0.520 

FPEO.  PPCO.  = 

0.0011440 

starting  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS: 

5 

22  9KORD 

13  84 

0.500  13 

84 

11 

81 

0.220 

11  61 

9 778 

0.100  9 

778 

5 

773 

0.400 

S 773 

1 494 

0.200 

FPEO.  PROO.  = 

0.0010560 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

22  OKOPD 

14  84  A 

• 

0.650  14 

84  A 

6 

773A 

0.090 

6 773A 

9 778 

0.210  9 

778 

1 

494 

0.100 

FREO.  PROD,  = 

0.0012285 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS: 

5 

22  BKORD 

14  84A 

0.650  14 

84A 

6 

773A 

0.090 

6 773A 

9 778 

0.210  9 

778 

3 

771 

0.200 

3 771 

1 494 

0.520 

FREO.  PROD.  * 

0.0012776 

START l\G  NODE 

IS: 

1 494 

NO. 

OF  pairs: 

5 

22  BKORO 

14  64  A 

0.650  14 

84  A 

6 

773A 

* 0.090 

6 773A 

9 778 

0.210  9 

770 

5 

773 

0.480 

5 773 

1 494 

0.200 

FREO.  PROD.  « 

0.0011794 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  pairs: 

4 

22  8<0RD 

14  04  A 

0.650  14 

04  A 

7 

774 

0.160 

7 77*. 

9 778 

• 

0.140  9 

778 

1 

494 

0.100 

FPEO.  PPOO.  * 

0.0014560 

STARTING  NOOE 

IS; 

1 494 

. NO. 

pf  PAIRS: 

S'  . . •' 

22  B«OPO 

14  84A- 

0.650  14 

04  A / * 

7 

774 

0.160 

7 774 

9 778 

0.140  9 

778 

3 

771 

0.200 

3 771 

1 494 

0.520 

• 

FREO.  PROD.  * 

0.0015142 

STARTING  node 

IS: 

1 494 

NO. 

OF  pairs: 

5 

22  B<ORD 

14  84  A 

0.650  14 

04  A 

7 

774 

0.160 

7 774 

9 778 

O.UO  9 

778 

5 

773 

0.480 

5 773 

1 494 

0.200 

FREO.  PROD.  * 

0.0013978 

starting  node 

IS: 

1 494 

NO. 

OF  pairs: 

5 

12  6K0R0 

14  84A 

0.650  * 14 

84A 

12 

81 A 

0.350 

12  8 1 A 

15  842 

0.190  15 

842 

3 

771 

0.060 

3 771 

1 494 

0.520 

FREO.  PROO.  * 

0.0013486 

STARTING  NODE 

IS:  ' 

1 494 

NO. 

OF  paips: 

4 

22  B<0«0 

14  84  A 

0.650  14 

04  A 

15 

842 

0.150 

15  842 

3 771 

0.060  # 3 

771 

1 

494 

0.520 

FREO.  PROD.  * 

0.0030420 
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TV* bio  3-6  (Continued) 


STARTING  NODE 

IS: 

I 494 

NO. 

OF  PAIRS: 

3 

22  6KOWD 

IS  842 

0.050 

15 

64*.  , 

3 

771 

0.060 

3 771 

1 494 

0.520 

FREQ.  P^OD.  * 

0.0015600 

USER  ENDING  NODE 

IS:.  23 

ICP 

STARTING  SOCE 

IS: 

l 494 

NO. 

OF  pairs: 

2 

?3  ICP 

11  81 

0.090 

11 

81 

1 

494 

0.080 

FREO.  PRCO.  a 

0.0072000 

• / 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  pairs: 

3 \ . V • 

23  ICP 

11  61 

0.090 

11 

81 

• 3 

771 

0.150 

3 771 

1 494 

0.S20 

. - # 

FREO.  PROD.  = 

0.0070200 

STARTING  NOOE 

IS: 

I 494 

NO. 

OF  PAIRS: 

3 

23  ICP 

11  81 

0.090 

11 

81 

5 

773 

0.070 

S 773 

♦ 1 494 

0.200 

FREO.  PROD.  a 

0.0012600 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

3 * 

23  ICP 

15  842 

0.070 

15 

842 

3 

771 

0.060 

3 771 

1 494 

0.520  * 

FREO.  PROD.  = 

0.0021840 

USER  ENDING 

NODE 

IS:.  24 

REFRO 

STARTING  NODE 

IS: 

• 1 494 

NO. 

OF  pairs: 

3 

24  REFRD 

2 722 

. 0.050 

2 

722 

3 

771 

0.070 

3 *771 

1 494 

0.520 

FREO.  PROD.  » 

0.0018200 

\ s • 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  pairs: 

2 

24  REFRD 

3 771 

0.080 

3 

771  . 

1 

494 

0.520 

FREO.  PROD.  = 

0.0416000 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

3 

24  REF  WO 

4 77  1 A 

0.140 

4 

77 1 A 

9 

778 

0.360 

9 776 

• 

1 494 

0.100 

FREO.  PROD.  = 

0.0050400 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4* 

24.  REFRD 

4 77 1 A 

0.140 

4 

77 1 A 

9 

778  * 

0.360 

9 778 

3 771 

0.200 

3 

771 

1 

494 

0.520 

FREO.  PROD.  = 

0.0052416 

starting  node 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

• 

24  REFRD 

4 77 1 A 

0.140 

4 

771 A 

9 

770 

0.360 

9 776 

5 773 

* 

0.480 

5 

773 

1 

494 

0.200 

FREO.  PROD,  a 

0.0048384 

STARTING  NOOE 

is: 

1 494 

NO. 

OF  pairs: 

4 

24  REFRD 

4 771 A 

0.140 

4 

77  1 A 

10 

778A 

0.460 

10  778A 

13  84 

0.320 

13 

84 

1 

494 

0.100 

FREO.  PROD,  a 

0.0020608 

STARTING  NOOE 

IS: 

1 494 

• NO, 

OF  PAIRS*. 

3..  . 

24  REFRD 

4 771 A 

0.140 

4 

77 1 A 

11 

81 

0.420 

11  61 

1 494 

0.080 

FREO.  PROO.  a 

0.0047040 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

24  REFRD 

4 77  1 A 

0.140 

4 

77 1 A 

11 

81 

0.420 

11  61 

3 771 

0.150 

3 

771  . • 

1 

494 

0.520 

FREO.  PROD.  = 

0.0045864 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS*. 

2 

24  REFRD 

5 773 

0.060 

5 

773 

1 

494 

0.200 

FREO.  PROD.  = 

0.0120000 

• 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS*. 

3 

24  REFRD 

6 773A 

0.060 

6 

773A 

9 

778 

0.210 

9 776 

1 494 

0.100 

FREQ.  PROD.  = 

0.0012600 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS*. 

4 

24  REFRD 

6 773A 

0.060 

6 

773A 

9 

770 

0.210 

9 778 

3 771 

0.200 

3 

771 

1 

494 

0.520 

FREO.  PROD.  * 

0.0013104 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

24  REFRD 

6 773A 

0.060 

6 

773A 

9 

778 

0.210 

9 778 

5 773 

0.480 

5 

773 

1 

494 

0.200 

FREO.  PROD.  = 

0.0012096 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS*. 

2 

24  REFRD 

13  84 

0.180 

13 

84 

1 

494 

0.100 

F»EO.  PROD.  = 

0.0180000 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  pairs: 

3 

24  REFRD 

13  84 

% 0.180 

13 

84 

5 

773 

0.090 

5 773 

1 494 

0.200 

FREO.  PROO.  = 

0.0032400 

. 

. 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS: 

3 

24  REFRD 

13  84 

0.180 

13 

84 

11 

81 

0.220 

11  61 

1 494 

0.080 

FPEQ.  PROD.  = 

0.0031680 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

24  REFRD 

13  64 

0.180 

13 

84 

11 

81 

0.220 

11  81 

3 771 

0.150 

3 

771 

1 

494 

0.520 
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Tabic  3-6  (Continued) 


FREQ.  PROD*  ^ 

0.0030868 

i 

starting  NODE 

IS? 

1 494 

NO.  OF  PAIRS: 

4 

t 

?4  R£FRD 

14  04  A 

0.350  14 

04  A 

15 

842 

0.150 

15  ‘6*2 

3 771 

. 0.060  3 

771 

1 

494 

0.520 

FREQ* * PROD*  = 

0.0016360 

• 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

o V • 

24  REFRO 

15  84? 

0.140  15 

84?  • 

3 

771 

0.060 

3 771 

1 494 

0.520 

. 

FREQ.  PROD.  - 

0.0043680 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

3 

24  REFRO 

17  872 

0.220  17 

872 

9 

778 

0.050 

9 778 

. 

1 494 

0.100 

FREQ.  PROD.  - 

0.0011000 

STARTING  NOOE 

IS: 

1 494 

NO.  OF  PAIRS: 

4 * 

i 

2u  R£FRQ 

17  872 

0.220  17 

872 

9 

778 

0.050 

9 778 

3 771 

0.200  3 

771 

1 

494 

0.520 

FR£0*  PROD.  * 

0.0011440 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

4 

* 

2 * REFRD 

17  072 

0.220  17 

07? 

9 

778 

0.050 

9 77b 

5 773 

0.480  5 

773 

1 

494 

0.200 

FREQ.  PROO.  = 

0.0010560 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

4 

Zu  REFRD 

18  872A 

0.550  18 

872A 

10 

778A 

0.060 

10  77oA 

13  0«. 

0.320  13 

84 

1 

494 

0.100 

FREQ.  PROD.  * 

0.0010560 

USER  ENOING 

► NOOE 

IS: . 25 

SPOTBUY 

% 

STARTING  NODE 

IS: 

I 494 

NO.  OF  PAIRS: 

3 

25  SPOTBUY 

2 72? 

0.180  2 

722 

3 

771 

0.070 

3 .771 

1 494 

0.520 

FREQ.  PROD.  = 

0.0065520 

* * 

starting  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

i * ' 

25  SP0T9UY 

2 722 

0.180  2 

722 

5 

773 

0.080 

5 773 

1 494 

0.200 

"■ 

FREQ.  PROD.  = 

0.0028800 

starting*  node 

IS: 

1 494 

NO.  OF  PAIRS*. 

2 

25  SP0T8UY 

3 771 

0.440  3 

771 

1 

494 

0.520 

FREQ.  PROD.  * 

0.2287999 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

4 . 

25  SPOTBUY 

4 77 1 A 

0.860  4 

7*71  A 

2 

722 

0.080 

2 722 

3 771 

0.070.  3 

771 

1 

494 

0.520 

FREQ.  PROD.  = 

0.0025043 

STARTING  node 

IS: 

1 494 

NO.  OF  PAIRS: 

4 

25  SP0T9UY 

4 7 7 1 A 

0.860  4 

77 1 A 

2 

722 

0.080 

2 722 

5 773 

0.080  5 

773 

1 

494 

0.200 

FREQ.  PROD.  = 

0.0011006 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

3 

25  SPOTBUY 

4 771 A 

0.e60  4 

771 A 

9 

778 

0.360 

9 778 

1 494 

0.100 

FREO.  PROD.  = 

0.0309600 

STARTING  NODE 

IS: 

1 494 

NO.  OF  pairs: 

25  SPOTBUY 

4 77 1 A 

0.660  4 

77 1 A 

9 

778 

0.360 

9 778 

3 771 

0.200  3 

771 

1 

494 

0.520 

FREO.  PROD.  « 

0.0321984 

STARTING  NODE 

IS,: 

1 494 

NO.  OF  PAIRS: 

4 

25  SPOTBUY 

4 7 7 1 A 

0.860  4 

771 A 

9 

776 

0.360 

9 778 

5 773 

0.480  5 

773  , . 

1 

494 

0.20Q 

FREO.  PROD.  = 

0.0297216 

SMARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

5 

. 25  SPOTBUY 

4 77 1 A 

0.660  4 

771 A 

10 

778A 

0.460 

10  778A 

6 773A 

0.480  6 

773A 

9 

778 

0.210 

9 778 

1 494 

0.100 

FREO.  PPOO.  = 

0.0039876 

starting  node 

IS: 

1 494 

NO.  OF  PAIRS: 

6 

25  SPOTBUY 

4 77 1 A 

0.860  4 

771 A 

10 

778A 

0.460 

10  770A 

6 773A 

0.480  6 

773A 

9 

778 

0.210 

9 778 

3 771 

0.200  3 

771 

1 

494 

0.520 

FREO.  PROO.  = 

0.0041471 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

6 

25  SPOTBUY 

4 77  1 A 

0.860  4 

77 1 A 

10 

778A 

0.460 

10  778A 

6 773A 

0.480  6 

773A 

9 

778 

0.210 

9 77  b 

5 773 

0.480  5 

773 

1 

494 

0.200 

FPEO.  PROO.  s 

0.0038201 

STARTING  NOOE 

IS: 

1 494 

NO.  OF  PAIRS: 

5 

25  SPOTBUY 

4 77  1 A 

0.860  4 

771 A 

10 

778A 

0.460 

10  778A 

6 773A 

0.480  6 

*77  3 A 

11 

61 

0.160 

11  81 

1 494 

0.080 

FREO.  PROD.  * 

0.0024306 

* 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS: 

6 

25  SPOTBUY 

4 771 A 

0.660  4 

771 A 

10 

770A 

0.460 

10  770 A 

6 773A 

0.480  6 

773A 

11 

81 

0.160 

11  81 

3 771 

0.1S0  3 

771 

1 

494 

0.520 

FREO.  PROO.  = 

0.0023690 

STARTING  NOOE 

IS: 

1 494 

NO.  OF  PAIRS: 

5 

• 

25  SPOTBUY 

4 771 A 

0.660  4 

77 1 A 

10 

778  A 

0.460 

10  .7-7  8 A 

7 774 

0.730  7 

774 

9 

778 

0.140 

9 778 

1 494 

0.100 
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Tr.ble  £-G  (Continued) 


FPEO.  PROO.  = 

» 

0.00*000 

> 

STARTING  NODE 

IS: 

1 49«* 

NO.  OF  PAIRS: 

6 

25  SROTtfUT 

* 7 7 1 A 

0.860  * 

77 1 A 

10 

776A 

0.460 

10  778 A 

7 77* 

0.730  7 

774 

9 

778 

0.140 

9 778 

3 771  . 

0.200  3 

771  V . 

1 

494 

0.520’ 

FREO.  PROO.  * 

0.00*20*8 

r 

STARTING  NODE 

IS: 

1 *9* 

NO.  OF  PAIRS: 

6* 

25  SPOT0UY 

* 77 1 A 

0.860  * 

771 A • 

10 

778A 

0.460 

10  778 A 

7 77* 

0.730  7 

77* 

9 

77B 

0.140 

9 778 

5 773  ‘ 

0 . *80  5 

773 

1 

494 

0.200 

FREO.  PROO.  * 

0.0038813 

+ 

STARTING  NODE 

IS: 

1 *9*. 

NO.  OF  PAIRS: 

4 

25  S°OTBUY 

* 77 1 A 

0.860  * 

771 A 

10 

778A 

0.460 

10  778A 

11  81 

0.110  11 

81 

1 

494 

0.080 

frc;q.  PROO.  * 

0.003*813 

starting  NODE 

IS: 

1 *9«* 

NO.  OF  PAIRS: 

5 

• 

25  SpOT9UY 

* 77 1 A 

0.860  * 

77 1 A 

10 

778  A 

0.460 

10  778A 

11  81 

0.110  11 

81 

3 

771 

0.150 

3 771 

1 *9* 

0.520 

FREO. • pROO.  * 

i 

0.00339*2 

starting  node 

IS: 

1 49* 

NO.  OF  PAIRS: 

4 

2S  SPOT9UY 

* 77  1 A 

0.660,  4 

77 1 A 

10 

778A 

0.460 

10  778 A 

13  8* 

0.320  13 

64 

1 

494 

0.100 

FREO.  PROO.  = 

0.0126592 

starting  node 

IS: 

1 *9* 

NO.  OF  PAIRS: 

5 ..  . 

25  SPOT0UY 

* 77 1 A 

0.660  * 

77 1 A 

10 

778A 

0.460 

10  7 78 A 

13  8* 

0.320  13 

84 

5 

773 

0.090 

5.  773 

1 *9* 

0.200 

FREO.  PROO.  2 

0.0022787 

STARTING  NODE 

IS:. 

1 *9* 

NO.  OF  PAIRS: 

5 

25  S°OT0UY 

* 77 1 A 

0.e60  4 

77 1 A 

10 

778A  ' 

0.460 

10  776A 

13  6* 

0.320  13 

84 

11 

81 

0.220 

11  81  * 

» • 

1 *9* 

0.080 

• 

FREC.  PROD.  = 

0.0022280 

i ; 

STARTING  NOOE 

IS! 

1 *9* 

NO.  OF  PAIRS: 

6 

25  SPOTBUY 

* 7 7 1 A 

0.860  4 

77 1 A 

10 

778A 

0.460 

10  7 76 A 

13  8* 

0.320  13 

84 

11 

61 

0.220 

11  61 

3 771 

0.150  3 

771 

1 

494 

0.520 

FREO.  PROD.  * 

0.0021723 

t * 

STARTING  NODE 

IS: 

1 *9* 

NO.  OF  PAIRS! 

5 • 

25  SPOTBUY 

* 771 A 

0.e60  4 

‘77 1 A 

10 

778A 

0.460 

10  7 78 A 

•15  8*2 

0.220  15 

8*2 

3 

771 

0.060 

3 771 

1 *9* 

0.520 

FREO.  PP 00.  * 

0.002715* 

• • 

f 

STARTING  NODE 

IS: 

1 *9* 

NO.  OF  PAIRS: 

3 

25  SPOTBUY 

* 77 1 A 

0.660  4 

771  A 

11 

81 

0.420 

11  B 1 

1 *9* 

0.080 

FREO.  PROO.  = 

0.0288960 

• , 

STARTING  NOOE 

IS: 

1 *9* 

NO.  OF  PAIRS: 

4 .*  1 

i 

25  SPOTBUY 

* 77 1 A 

0.660  4 

‘77 1 A 

11 

81 

0.420 

11  81* 

3 771 

0.150  3 

771 

1 

494 

0.520 

FREO.  PROD.  = 

0.0281736 

STARTING  NOOE 

IS: 

1 *9* 

NO.  OF  PAIRS: 

4 

25  SPOTBUY 

* 77 1 A 

0.660  4 

77 1 A 

11 

81 

0.420 

11  81 

5 773 

0.070  5 

773 

1 

494 

0.200 

FREO.  PROO.  » 

0.0050568 

• 

starting  node 

IS: 

1 *9* 

NO.  OF  PAIRS: 

4 

25  SPOTBUY 

* 77 1 A 

0.860  4 

771 A 

11 

61 

0.420 

11  61 

9 778 

0*100  9 

778 

‘ 1 

494 

0.100 

FREO.  PROO.  » 

• 

0.0036120 

STARTING  NODE 

IS: 

1 *9* 

NO.  OF  PAIRS: 

5 

25  SPOTBUY 

* 77 1 A 

0.860  4 

771 A 

11 

81 

0.420 

‘ 11  81 

9 778 

0.100  9 

778 

3 

771 

0.200 

3 771 

1 *9* 

0.520 

* • 

FREO.  PROO.  = 

0.0037565 

STARTING  node 

IS: 

1 *9* 

NO.  OF  PAIRS: 

5 

25.  SPOT9UY 

* 77 1 A 

0.660  4 

’ 77 1 A • ' 

*11 

81 

0.420 

11  81  # 

9 779 

0.100  9 

N 776 

5 

773 

' 0.480 

5 773 

1 *9* 

0.200 

\ • 

FREO.  PROO.  * 

0.003*675 

STARTING  NOOE 

IS: 

1 *9* 

NO.  OF  pairs: 

5* 

25  SPOTBUY 

* 77 1 A 

0.860  4 

77 1 A 

12 

81 A 

0.510 

12  6 1 A 

10  778A 

0.130  10 

77eA  - 

13 

84 

0.320 

13  6* 

1 *9* 

0.100 

FREO.  PRCO.  * 

0.00182*6 

; 

• 

STARTING  NODE 

IS: 

1 49*. 

NO.  OF  pairs: 

5* 

25  SPOTBUY 

* 77 1 A 

0.860  4 

77 1 A 

12 

au 

0.510 

12  8 1 A 

15  8*2 

0.190  15 

842 

3 

771 

0.060 

3 771  1 494  * 0.520 

FREO.  PROO.  ■ 0.0026000 
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Table  3-G  (Continued) 


starting  node 

IS  t 

1 494 

NO. 

OF  pairs: 

4 

25  SPOTBUY 

4 7 7 1 A 

o.ebo 

4 

77 1 A 

15 

842 

0.060 

15  6<.2 

3 771 

0.060 

3 

771 

i 

494 

0.520 

FPEO.  PPOO.  a 

0.0016099 

i 

STARTING  NOOE 

1ST 

1 494 

NO. 

of  pairs: 

2 

25  SPOTBUY 

5 773 

0.220 

5 

773 

i 

494 

0.200 

FPEO.  PROD.  c 

0.0440000 

..  , 

STARTING  NODE 

IS« 

1 494 

NO. 

OF  pairs: 

3 

25  SPOT3UY 

6 773A 

0.220 

6 

773A 

9 

776 

0.210 

9 778 

1 494 

0.100 

FRED.  PP 00.  * 

0.0046200 

starting  node 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

25  SPOTflUY 

6 773  A 

0.220 

6 

773A 

9 

778 

P 

0.210 

9 776 

3 771 

0.200 

3 

771 

1 

494 

0.520 

FREO.  PPCO.  » 

0.0048048 

• 

5TARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

25  SP0T8UY 

6 773A 

0.220 

6 

773A 

9 

778 

0.210 

9 776 

5 773 

0.480 

5 

773 

1 

494 

0.200 

FREQ*  P POO.  a 

0.0044352 

STARTING  node 

IS: 

1 494 

NO. 

OF  pairs: 

3 

25  SP0T6UY 

6 773A 

0.220 

6 

773A 

11 

81 

0.160 

11  81 

1 494 

0.080 

FPEQ.  PROO.  = 

0.0028160 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 

25  SPOTBUY 

6 773A 

0.220 

6 

773A 

* 11 

81 

0.160 

11  81 

3 771 

0.150 

3 

771 

1 

494 

0.520 

FPEO.  PROD.  = 

0.0027456 

STARTING  NOOE 

IS: 

1 494 

NO, 

OF  PAIRS: 

3 

25  SPOT0UY 

7 774 

0.110 

7 

774 

9 

778 

0.140 

9 778 

1 494 

0.100 

FPEQ.  PROD.  = 

0.0015400 

* 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  pairs: 

4 

i 

25  SPOTBUY 

7 774 

0.110 

7 

774 

9 

770 

0.140 

9 778 

3 771 

*0.200 

3 

771 

1 

494 

0.520 

FPEO.  P POO.  * 

0.0016016 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

4 • 

25  SPOTBUY 

7 774 

0.110 

7 

774  • 

9 

778 

0.140 

9 778 

5 773 

0.480 

5 

773 

1 

494 

0.200 

FREQ.  PPOO.  » 

0.0014794 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

6 

25  SPOTBUY 

8 774A 

0.550 

6 

774A 

10 

770A 

0.180 

10  778A 

6 773A 

0.480 

6 

773A 

9 

776 

0.210 

9 778 

3 771 

0.200 

3 

771 

1 

494 

0.520 

FPEO.  PROD.  * 

0.0010378 

■ 

STARTING  NODE 

1ST 

1 '494 

too. 

OF  PAIRS: 

5 

25  SPOTBUY 

8 7 74  A 

0.550 

6 

774A 

lb 

778A 

0.180 

10  *776  A 

7 774 

0.730 

7 

774 

9 

778 

0.140 

9 . 778 

* 

1 494 

• 

0.100 

’» 

FREO.  PPOO.  « 

0.0010118 

* i • * 

STARTING  NOOE 

is: 

1 494 

NO, 

OF  PAIRS: 

v6  / * 

25  SPOTBUY 

8 774A 

0.550 

8 

774A 

10 

770A 

0.180 

10  778A 

7 774 

0.730 

7 

774  • 

9 

778 

0.140 

9 778 

3 771 

0.200 

3 

771 

1 

494 

0.520 

FREQ.  PROD.  * 

0.0010522 

• 

STARTING  NOOE 

1ST 

1 4Q4 

NO. 

OF  PAIRS: 

• 

4 " 

25  SPOTBUY 

8 774A 

0.550 

6 

774A 

10 

778A 

0.160 

10  7 7b A 

13  84 

0.320 

13 

84 

1 

494 

0.100 

FPEO.  PROD.  = 

0.0031O80 

» 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS: 

2 

25  SPOTBUY 

9 778 

0.090 

9 

778 

1 

494 

0.100 

rREQ.  PROD.  * 

0.0090000 

STARTING  NODE 

IS: 

1 494 

NO. 

OF  PAIRS: 

3 

25  SPOTBUY 

9 778 

0.090 

9 

778 

3 

771 

0.200 

3 771 

1 494 

0.520 

FREO.  PROO.  * 

0.0093600 

STARTING  NOOE 

IS: 

1 494 

NO. 

OF  PAIRS? 

3 

25  SPOTBUY 

9 778 

0.090 

9 

778 

5 

773 

0.480 

5 773 

1 494 

0.200 

*RE3.  PPOO.  a 

0.0006400 

„ , 

STARTING  NOOE 

IS: 

1 494 

'no. 

OF  PAIPS: 

4 

25  SPOTBUY 

10  778A 

0.110 

10 

77eA 

6 

773A 

0.480 

6 773A 

9 778 

0.210 

9 

778 

1 

494 

0.100 

rREO.  PROO.  a 

0.0011088 

. 

STARTING  NOOt 

IS: 

1 494 

NO. 

OF  PAIRS: 

5 

25  SPOTBUY 

10  7 78A 

0.110 

10 

778A 

6 

773A 

0.400 

6 773A 

9 778 

0.210 

9 

778  * 

3 

771 

0.260 

3 771 

1 494 

0.520 

• 

FREO.  PPOO.  « 

0.0011532 

* 
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Tabic  3-G  (Continued) 


STARTING  node 

IS: 

1 494 

NO.  OF  PAIRS:  5 

. 

25  SPOTBUY 

10  778A 

0.110  10  77A A 

6 

773A 

0.480 

6 773A 

9 778 

0.210  9 778 

5 

773 

0.4B0 

5 773 

1 494 

0.200 

rREO.  PROD,  = 

0.0010644 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS:  4 

25  SPOTBUY 

10  778A 

0.110  10  778A 

7 

774 

0.730 

7 774 

* 

9 778 

0.  140  9 778 

1 

494 

0.100 

FREQ.  PROD*  = 

0.0011242 

r 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS:  5* 

25  SPOTBUY 

10  7 78A 

0.110  10  778A 

7 

774 

0.730 

7 774 

9 778 

0.140  9 778 

3 

771 

0.200 

3 771 

1 494 

0.520 

FREO.  PROD • s 

0.0011692 

STARTING  NOOE 

IS: 

1 494 

NO.  OF  PAIRS:  5 

25  SPOTBUY 

10  778A 

0.110  10  778A 

7 

774 

0.730 

7 774 

9 778 

0.140  9 778 

5 

773 

0.480 

5 773 

1 494 

* 0.200 

FREO.  PROD.  = 

0.0010792 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS*.  3 

25  SPOTBUY 

10  778 A 

0.110  10  778 A 

13 

84 

0.320 

13  84 

1 494 

0.100 

FREO.  PROD.  = 

0,0035200 

STARTING  NODE 

IS: 

1 494 

NO.  OF  pairs:  4 

25  SPOTBUY 

20  873A 

0.530  20  873A 

10 

778A 

0.060 

10  778 A 

13’  84 

0.320  13  84 

1 

494 

0.100 

FREO,  PROD.  «= 

0.0010176 

USER  ENDING 

• NODE 

IS: ♦ 26 

CANC 

STARTING  NODE 

IS: 

1 494 

NO.  OF  PAIRS:  3 

26  CANC 
3 771 

2 722 

1 494 

0.070  2 722, 

0.520 

3 

771 

0.070 

FREO.  PROD.  = 

0.0025480 

STARTING  NODE 

IS: 

1 494  ' 

NO.  OF  PAIRS:  3 

26  CANC 
5 773 

? 722 

1 494  ♦ 

0.070  2 . 722 

0.200  ^ 

5 

773 

0.080 

FREO.  PROD.  = 

0.0011200 

\ .•  , 

• 
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CONCLUDING  REMARKS 

The  purpose  of  the  research  effort  reported  on  here  was  to  apply  a methodology 
which  would  document  the  relationship  between  requisition  processing  throughput 
times  and  the  resources  dedicated  to  requisition  processing.  The  major  result  of 
the  study  is  a reinforcement  of  an  insight  gained  earlier  in  the  development  of  the 
method  (and  concluded  by  other  studies): 

that  requisition  processing  throughput  times  are  sensitive  to  many  factors 
other  than  simply  manpower  available  to  do  the  processing. 

As  a result,  available  data  reflect  so  many  influences  on  elapsed  processing 
times  that  those  portions  of  throughput  time  distributions  which  are  attributable 
to  available  man-hours  are  difficult  to  document.  The  design  of  a statistical  model 
to  estimate  the  influence  of  resource  constraints  would  presumably  require  a sub- 
stantial amount  of  longitudinal  performance  data  as  well  as  data  describing  the  per- 
haps complex  priority  structure  governing  the  allocation  of  resources  between 
requisition  processing  and  other  tasks.  It  is  not  clear  that  existing  data  sources 
would  support  such  a statistical  model. 

Further,  it  is  not  clear  that,  practically  speaking,  requisition  processing  as 
currently  practiced  is  in  fact  resource  constrained  at  all.  It  can  be  argued  that  a 
legitimate  and  necessary  task  of  the  Navy  Supply  System  is  the  maintenance  of  ex- 
cess capacity  against  potential  emergency  requirements.  If  this  is  the  case,  then 
current  performance  data  will  fail  to  reveal  a sensitivity  of  throughput  time  dis- 
tributions to  changes  in  the  work  load/capacity  relationship.  The  study  of  resources 
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assigned  to  processing  and  billing  requisitions  might  appropriately  be  directed  to 
the  size  and  expense  of  excess  capacity.  Additionally,  structural  issues  might 
arise.  If  excess  capacity  is  indeed  being  maintained,  it  would  appear  to  be  un- 
evenly maintained.  Requisition  processing  resources  appear  to  be  under-utilized 
to  a greater  degree  than  materials  handling  resources.  In  fact  the  issue  of  the 
supply  system's  maintaining  a capability  might  lead  to  performance  concepts  some- 
what different  from  those  imposed  on  the  system  from  the  standpoint  of  its  proc- 
essing of  requisitions. 

Hie  results  of  this  research  provide  a model  of  requisition  processing  which 
is  not  firmly  verified  from  observed  requisition  processing  practices . Ignoring, 
generically  speaking,  the  tasks  and  the  training  in  tasks  which  comprise  requisi- 
tion processing,  it  does  not  appear  that  the  distributions  of  throughput  times  associ- 
ated with  requisition  processing  arise  in  general  as  the  result  of  resource  constraints. 
However,  the  ability  to  parameterize  the  model  provides  great  flexibility  in  repli- 
cating those  subportions  of  processing  practices  which  might  be  determined  to  be 
resource-constrained,  or  in  extrapolating  the  implications  of  work  load/capacity  re- 
lationships far  beyond  circumstances  currently  observed. 
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APPENDIX  A 


• multipriority  queueing  model 

• Kolmogorov- Smirnov  test 

• results  for  nonpreemptive  queues 


A-l 


This  puper  pri.  • <■.  : llu  mathei.uitieal  models  used  to  study  the  queueing  theoretic 
acprrp;  of  the  throughput,  lime  :t.«so<*i-led  with  requisition  processing  and  materials 
handling'.  A model  was  prepared  for  circumstances  of  each  of  one,  two,  or  three 
priorities.  The  underlying  input  stream  for  each  priority  is  assumed  to  be  Poisson 
in  all  of  the  models.  This  is  clearly  not  an  exact  representation  of  reality,  but  it 
can  be  shown  that  the  more  realistic  assumption  that  inputs  arc  at  constant  inter- 
vals in  batches  whoso  random  sizes  arc  Poisson  can  be  well-approximated  by  a 
Poisson  process.  This  is  demonstrated  through  the  following  argument: 

Pr  | a total  of  n units  have  arrived  on  or  before  the  end  of  the  l;-’1  uniformly 
spaced  interval} 

---  Pr  | 1;  modules  of  random  and  Poisson-distributed  size  lead  to  a total 
ol  11} 

- Pr  | sum  of  1;  Poissons  totals  n } . 


But  it  is  well  known  that  1:  identical  Poissons  with  parameter  a sum  to  another  Pois- 
son, Il'r>  one  with  parameter  Hence, 

Pr  j sum  ofl;  Poissons  totals  n} 

-la>  j} 

~ c (ha)  /n!, 


that  is,  Poisson  with  mean  lea.  In  view  of  the  facts  that  the  system  will  not  be  able 
to  react  immediately  to  the  input  and  that  the  time  between  batches  can  be  lowered, 
the  overall  process  (not  just  that  evaluated  at  a point  of  arrival)  will  largely  act  as 
a Poisson  process  even  though  the  proved  result  is  not  valid  for  times  between 
batches. 

A routine  has,  in  fact,  been  provided  (under  the  name  of  POIS)  for  testing  whether 
any  particular  set  of  data  are  Poisson  distributed.  In  addition,  since  the  models  also 
assume  that  service  times  are  exponentially  distributed,  a second  routine  (names  EXP) 
has  been  written  for  the  testing  of  exponentiality. 

Generally,  the  easiest  and  most  familiar  way  to  test  for  Poisson  or  exponential 
2 

character  is  to  use  a X goodness -of-f it  test  on  the  data  presented  in  block  histogram 
form  against  a theoretical  distribution  with  each  parameter  replaced  by  its  maximum- 
likelihood  estimator,  which  is 
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n ^ - At 

X = n/  £ t.  for  both  the  exponential  Ae  and  Poisson  (At)ne  '/n  ! , 

i = 1 1 •' 

where  t.  is  the  time  between  the  (i  - l)St  and  i^1  occurrences.  The  resulting  sta- 
1 l 

tistic  is  then 


vi  I'wX1' 

1 = 1 

th 

where  o.  is  the  number  observed  in  the  i frequency  class  (out  of  a total  of  between 

. th 

10  and  20  classes),  e.  the  number  expected  in  the  i frequency  class  if  the  hypothe- 
sized distribution  were  correct,  and  k the  number  of  degrees  of  freedom,  always 
equal  to  the  total  number  of  classes  less  one  and  then  minus  one  for  each  parameter 
estimated.  Of  course,  the  usual  px-ecautions  must  be  taken  to  keep  the  number  in 

any  class  from  being  too  small  (a  rule  of  thumb  being  less  than  five). 

2 

Great  care  should  always  be  exercised  in  doing  X goodness -of-fit  tests  and  the 

analyst  would,  of  course,  be  well  advised  to  search  for  a definitive  exposition  on  the 

2 

subject  in  the  statistical  literatui’e.  The  basic  weaknesses  of  the  X test  are  its  re- 
quii'ement  for  large  samples,  its  heavy  dependence  on  the  choice  of  the  number 
and  position  of  the  time-axis  intervals,  and  its  possibly  very  high  Type  II  error  (this 
is  expressed  in  terms  of  the  probability  of  accepting  a false  hypothesis)  for  feasible 
alternative  distributions.  In  view  of  these  difficulties,  we  are  instead  going  to  sug- 
gest two  tests  for  use  in  our  context:  the  Kolmogorov-Smimov  (K-S)  for  Poisson 
fits  and  the  F-test  for  exponentials. 

The  K-S  test  compares  deviations  of  the  empirical  CDF  from  die  theoretical 
CDF,  and  uses  as  its  test  statistic  the  maximum  absolute  deviation,  that  is. 


E = maximum 


.th 


n.  - F (N.)  , 

j y I 


where  m is  the  j ordered  (ascending)  observation,  and  F (N^.)  is  the  Poisson  CDF, 


1. 


critical  values  can  be  found  in  almost  any  statistics  text. 
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Z [e"A  A1/!!],  with  A.  = 1/t , t the  sample  mean.  The  5 percent  and  1 percent 
i = 0 

critical  values  for  the  K-S  are  stored  in  the  computer,  and  the  hypothesis  of  "Pois- 
sonness"  is  rejected  then  if  the  value  of  E exceeds  the  tabulated  critical  value. 

Statistics  are  given  in  Table  1. 

This  test  was  originally  derived  for  fitting  continuous  CDFs,  but  can  be  used  as 
a slightly  more  conservative  test  in  discrete  cases,  with  its  power  improving  with 
increasing  sample  sizes.  ( ..  , 

For  the  exponential  F test,  the  first  r (=  n/2)  and  then  (n  - r)  of  a set  of 

n hypothesized  exponential  interoccurrence  times  |t  1 arc  grouped,  and  S.  is  used 

th  1 1 

to  denote  the  i normalized  spacing,  that  is, 

S.  = (n  = i+  1)  (t.-  <tt-  1))  (tQ3  0). 

Then  the  j S. } are  independent  and  identically  distributed  exponentials  with  exactly 
the  same  mean  as  the  underlying  distribution.  Thus  it  follows  that  the  quantity 

£ S,/V 

F = _i- 

£ S,./(n  - r) 

i = r + 1 

is  the  ratio  of  two  gammas  and  is  distributed  as  an  F distribution  with  2 r and  2 (n  - r) 
degrees  of  freedom  where  the  hypothesis  of  exponentiality  is  true.  Therefore,  a two- 
tailed  F test  would  be  performed  on  the  F statistic  calculated  from  the  data  in  order 
to  determine  whether  the  stream  is  indeed  truly  exponential.  The  left  and  right  F 
critical  points  for  a and  b degrees  of  freedom  at  the  5 percent  level  of  significance 
(say,  F _ (a,  b)  and  F (a,  b),  respectively)  can  be  found  from  the  following 
approximate  formula  (tested  to  be  within  0. 6 percent  accuracy  of  exact  values): 


A -4 


b - 3.986 
. 1414  b - . 2864 


_ . , . 3 t 1«  739 

.975  (a’  .1197  a + .1108 


< 


+ . 145  - . 00170  a + 


.06150  a - 2.706 
b+  30 


^ F.025  (a’b)  = l/F.975  (b’a)*1 


This  approximate  approach  is  a great  help  in  computer  applications  since  the  storage 
of  a complete  F table  can  be  replaced  by  the  use  of  these  equations  and,  further,  no 
interpolation  formula  is  needed  as  might  be  required  by  the  F-table-storage  method. 


1.  See  Table  1. 
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Table  1 


KOLMOGOROV-SMIRNOV  CRITICAL  VALUES 


Sample  Size 

5 Percent 

1 Percent 

5 

.56 

.67 

6 

.53 

.63 

7 

.50 

.60 

8 

.47 

.56 

9 

.44 

.53 

10 

.41 

.49 

11 

.40 

.47 

12 

.38 

.45 

f 13 

.37 

.43 

I 14 

.35 

.41 

15 

.34 

.40 

16 

.33 

.39 

17 

.32 

.38 

18 

.31 

.37 

19 

.30 

.36 

20 

.29 

.36 

21 

.29 

.35 

22 

.29 

.34 

23 

.28 

.33 

24 

.28 

.32 

25 

.27 

.32 

26 

.27 

.31 

27 

.26 

.31 

28 

.26 

.30 

29 

.25 

.30 

H 30 

.25 

.29 

^30 

1.36/Vn 

1. 63/Vn" 
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So,  for  the  marginal  queueing  model  implemented  at  each  channel,  if  the 
Kolmogorov-Smimov  and  F statistics  do  not  lead  to  rejection,  it  is  assumed  that  work 
units  arrive  as  a Poisson  process  to  a single  exponential  channel  and  that  upon  arrival 
to  the  system  each  unit  is  designated  to  be  a member  of  one  of  three  priority  classes 
(or  less  as  circumstances  dictate).  The  usual  convention  is  to  number  the  priority 
classes  so  that  the  smaller  the  number,  the  higher  the  priority.  Let  it  further  be 
assumed  that  the  arrivals  of  the  first  or  highest  priority  have  mean  arrival  rate  of  a^ 
work  units  per  unit  time,  that  the  second  or  middle  priority  units  have  mean  rate  a0 
work  units  per  unit  time,  and  that  the  third  or  lowest  priority  units  have  mean  a work 

u 

units  per  unit  time,  such  that  their  sum  is  called  a.  The  corresponding  service  x-ates 
shall  then  be  u , u 2,  and  uy  work  units  per  time  unit  for  prioi-itics  1,  2,  and  3,  re- 
spectively. Let  it  further  be  supposed  that  the  fii'st  priority  items  have  the  right  to 
be  sei’vcd  ahead  of  the  others,  but  that  once  a service  of  a priority  2 or  3 work  unit 
is  begun,  it  cannot  be  interrupted  by  preemption. 

In  light  of  these  assumptions,  it  has  been  shown  [see  Cobham  (1954)  or  Morse 
(1958)]  that  the  expected  number  of  work  units  in  the  queueing  system  for  each  priority 

can  be  faii-ly  easily  found  in  terms  of  the  input  and  service  parametci's.  If  Q , Q , 

X z 

and  are  used  to  denote  tht.se  averages,  then  urc  have  [see  Equation  (A. 4)  in  the 
Appendix] : 


3 


3 


k = 1 


a 


2 


Q, 


2 (1  - a1/u]L)  (1  - a^Uj^  - a2/u2)  + 


(1) 


u 


2 


3 


Q 


k - 1 
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The  mean  system  waiting  times,  say  W(l),  W(2),  and  W(3),  are  then  found  by 
applying  Little's  formula,  Q - aW,  on  Equation  (1),  so  that  \V(1)  = Q^/a  , W(2) 

= Q /a  , and  W(3)  = Q /a  . The  total  average  system  wait  can  then  be  obtained 

A Z ti  O 

by  the  weighted  average  of  \V(1),  W(2),  and  \V(3),  namely, 


The  variances  of  the  system  delays  for  the  three  priorities  can  be  determined 
using  results  of  Kesten  and  RunneaLurg  (1037)  which  were  derived  in  a manner 
very  similar  to  the  work  of  Cobham.  Without  going  into  the  details  of  the  required 
results,  suffice  it  t«  say  that  the  key  ones  are  those  which  give  the  second  moments 
of  the  waits  in  line  as 


W = (a1/a)  W(l)  + (a2/a)  \V(2)  + (a3/a)  W(3). 


3 


3 


2 


(1  “ aj/Uj)' 


3 


+ 


+ 
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To  be  move  exact,  it  is  assumed  Hint  the  right-hand  inequality  in  (2)  is  binding 


and  thus  that 

Pr  ) | X - K[X]  | ^ k a j = l/\? . 

Now  assuming  further  that  the  probability  distribution  has  equal  probability  on 
cacli  side  of  the  mean, 

Pr  { X - E[X]  n k <r  ( =1/(2  l:2) 

and 

2 1 

Pr  { K[X]  -X>kir|  = 1/(2  k ). 

2 

So,  given  the  moan  E[XJ  and  the  variance  c (or,  equivalently,  the  standard  devi- 
ation a),  the  distribution  function  may  be  reconstructed  by  varying  k in  reasonably 
small  step  over  an  appropriate  range.  In  the  program  written  for  the  analysis 
this  is  done  automatically  for  each  subsystem  and  then,  for  any  specific  values  of 
the  input  parameters,  summary  information  about  the  queue  is  printed  out  in  the  form 
of  the  average  number  of  units  or  each  priority  in  the  system,  the  total  average  sys- 
tem wait,  the  variance  of  the  system  wait  for  each  priority,  the  (approxin.  .e)  proba- 
bility distribution  for  the  three  system  delays. 

The  distributions  for  the  thn.  • priorities  must  then  be  combined  in  order  to 

obtain  ’ e probabilities  for  the  to.nl  process.  This  is  done  by  the  usual  mixing  pro- 

th 

ceduro  s follows.  If  the  individual  prob;  tities  for  the  k priority  are  denoted  by 
j Pj  (k)»  3 ^ i 20 1 , and  the  combined  distribution  by  |C.,  l;f  i < 40  } , then 

C.  = (ax/a)  p.(l)  + (a2/a)  p.(2)  + (a3/a)  p.(3). 

For  circumstances  imder  which  a two-prioritj'  or  a no»pviority  service  protocol 
is  observed,  appropriately  amended  programs  have  been  provided  NAVSUP.  In  the 
following  appendix  may  be  found  additional  technical  discussions  of  queues  with  many 
priori!  s. 


1.  This  assumption  may  be  changed  by  altering  the  value  of  an  appropriate 
parameter  in  the  program. 

4 


« 
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RESULTS  FOR  NON  PREEMPTIVE  EXPONENTIAL  QUEUES 


WITH  MANY  PRIORITIES 


RESULTS  FOR  N ON PRE E IvIPTIV E EXPONENTIAL  QUEUES 
WITH  MANY  PRIORITIES 

As  was  noted  earlier  in  the  previous  discussion,  the  determination  of  the  mean 
system  sizes  and  waiting  times  can  be  found  via  a direct  expected-value  procedure 
such  as  that  due  to  Cobhatn  [1954]  or  through  the  more  classical  differential-differencc- 
equation  method  such  as  that  found  in  Morse  [195S].  Cobham's  approach  to  the  mean 
delay  for  display  was  selected  since  it  is  more  straightforward  and  can,  in  fact,  be 
just  as  well  used  to  handle  multiple  priorities  as  two  priorities. 

To  begin,  suppose  that  items  of  a priority  (the  smaller  the  number,  the 
higher  the  priority)  arrive  before  a single  channel  according  to  a Poisson  distribu- 
tion with  parameter  a^  (k  = 1, 2, . . . , p)  and  that  these  work  units  wait  on  a first-come, 
first-served  basis  within  their  respective  priorities.  Let  the  service  distribution 
for  the  k^1  priority  be  exponential  with  mean  l/u  . Whatever  the  priority  of  an  <t 

u 

in  sendee,  it  completes  its  service  before  another  item  is  admitted. 

Begin  by  defining 

rk = V\  (1 2 1:2  p> 

and 

(R<f°-Bp3R)- 

The  system  is  stationary  for  R^  = R<1. 

Then  suppose  that  a work  unit  of  priority  i arrives  at  the  system  at  time  t^  and 
enters  service  at  time  t, . Its  lim.  wait  is  thus  T = t - t_.  At  t.  assume  that  there 
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arc  n work  units  of  priority  1 in  the  line  ahead  of  this  new  arrival,  n of  priority  2, 

X 2 

of  priority  3,  etc.  Let  be  the  time  required  to  finish  the  item  already  in  serv- 
ice and  Sj_  be  the  total  time  required  to  serve  nj..  During  the  new  woi'k  unit's  wait- 
ing time  T , (say)  n'  items  of  priority  k < i will  arrive  and  go  to  service  ahead  of 

q ^ ^ 

this  current  arrival.  If  is  the  total  service  time  of  all  the  n^,  then  it  can  be  seen 
that 

V.E.  S,'(+,E.  w 


k = 1 


k = 1 


If  expected  values  are  taken  on  both  sides  of  the  foregoing,  then  we  find  that 

i - 1 i 


Wq  “ ' ),  Et\J  + E,  EtSkJ  + E'S0>- 

K **  1 h.  — JL 


Since  R.  , -=  R.  for  all  i,  R <1  implies  that  R,  , «=  1 for  all  i. 
i - 1 i i - 1 

To  find  E[Sq]  , observe  that  the  combined  service  distribution  is  the  mixed 
exponential,  which  is  formed  from  the  law  of  total  probability  as 


where 

■-£v 

The  random  variable  "remaining  time  of  service, " S^,  has  the  value  0 if  the  system 
is  idle  and  hcuce 

E[Sq]  = Pr  {system  is  busy}  E[SQ  Jbusy  systeml. 

But  the  probability  that  the  sj'stcm  is  busy  is 


p -V 

B(t)=.]T  ak(1"e  )/a 


k - 1 
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a • (expected  service  It'  ')  = a E (nk/a)  (l/iy 


= R, 


and 


E[Sq  j system  busy] 


= E[S0  | system  busy  with  k type  work  unit].  Pr  { work  unit  has  priority  k } 
Pi 

■ <1/uk>  <rk/K)- 

Therefore 


E[S0]  = 11  Yj  (lAlk}  (VR) 

k*=  1 


= X,  < w- 

k = 1 


(A.l) 


Since  n^  and  the  service  times  of  individual  work  units,  S^,  are  independent, 
E[S01  E(Dk  S«, 

• “ E!nkl  <1 


■ E[nkl/Uk- 


Utili.  ing  Little's  form,  la  then  gives 

E[S.  ] = a.  W^Vu. 

L k He  q k 


= r,  W^. 

k q 
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• Similarly, 


E[s'l  - i-Wl  ,, 

K K iv 


and  then  utilizing  the  uniform  property  of  the  Poisson  we  have 


Ei\i  - \ % Vuk- 


Therefore 


i - 1 


VV 


» . W<‘>  Y'  r + Y'  r + E[S  ] , 

11  q k^l  k k=x  k » 0‘ 


or 


£.  rk  Wq°  * E‘S0> 


w(0  „ !i_i 

q 1 - R 


i - 1 


(A.  2) 


The  solution  to  Equation  (A.  2)  was  found  by  Cobham,  after  whom  much  of  this  analysis 
follows,  by  induction  on  i,  after  a general  pattern  c:  merged  upon  iteration.  That 
solution  is 


w<‘>  = 


E[S01 


q (1  - Ii.  _ j)  (1  - R.) 


Using  Eolation  (A.  1)  finally  gives 


(i)  k = 1 
Wv ) = 


E.  <Vv 


q (1-R.^)  (1-R.) 


(A.  3) 


Note  that  (A.  3)  holds  as  long  as  R = r 1.  Of  couit»c,  the  individual  mean  sys- 

k = 1 k 

tem  delay  for  priority  i is  therefore 

= \V^.  + l/u.  • 

Q i 
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Therefore,  from  Little's  formula,  t’.o  m an  number  of  work  units  of  priority  i pres 
ent  in  the  system  is  given  by 


Q.  = a.  W 
1 i q 


,'(i>  ♦ 


a./u. 

i 1 


P 


ai  X (ri/v 

k = 1 

(1  - It.  _ x)  (1  - It.) 


and  that  the  total  expected  system  size  is 


(A.  4) 


Expressions  very  similar  to  that  of  Equation  (A.  3)  were  found  for  the  higher 
moments  of  the  line  delays  for  each  priority  by  Kesten  and  Runneuburg  (1957).  These 
results  can  also  be  found  in  the  more  readily  available  reference  Cohen.  (19G9).  The 
formulas  are  a bit  lc  gthy  and  will  not  be  directly  noted  here  but  instead  may  be  found 
within  the  program  in  the  calculations  leading  to  the  system  waiting  time  variance. 
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APPENDIX  B 
PROGRAM  LISTINGS 

• Multipriority  Queueing  Model 

• Kolmogorov^ Smirnov  Test  Statistic  for  Poisson  Fit 

• General  Kolmogorov- Smirnov  Test  Statistic  Routine 

• Convolution  Routine 

• Exponential  Fit  Test 
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MUL TTPRIORIT  Y QUEUEING  MODEL 
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1 >.  ihlb  1 ij  KhLChAf".  LLc.i'r-1  FLh  ihi.  3 I IGnl'lY,  3 A«u>lvf'.  I.'-U, 
1»  3 bLhVILL  KAIL  KhLELEK 

3 Lie-  AC'  1 ) , 1 C3,  M >,KC3,  M >,i\  C 3,  M >,  M 3,  a 1 ) , 2Cb  1 ) , 1 C 3,  a 1 ) 

a Llh  UC3,  ' 1 ),  VC  3),V.  (3),  AC  3)  , L<  3^  l ) 

5 L IK  L < 3 , *•  1 ) 

b * u 1 , l2,  3 L 3 AKL  IKK  bLhVlCL  KAl  ts»  L.lAEb  IV,  20, 'Ll- 
1 b*  IF  1K.LKL  Af%L  t\C  Kh.iLhl'lltb  YbUbhCLLL  L'bE  KhCCnttiv.  Lu?Ci\E 
lb  u a=  1 • 30b 

19  L 1 = l . 1 * LA 
LG  LL  = 2*Lil 
LI  L 3 = L 2 

22*  M,A2,'.A3  AKL  1 KL  AKhlVAL  KA'lbb,  Lli\Lb  23,2', ''Lb 
L 3 M = .3Vl 
LA  AL= • bbL 
Lb  A 3=.7LV 

L6  h 1 = A 1 / U 1 + AL/U2  + A3/U3 
L 7 IF  h 1 > 1 • CL  1L  LVu 
Lb  A=A1+AL+A3 
LV  31  = 1 .-A  1/Ll 
3G  LL=M-AL/LL 

31  b3=bL-A3/L3 

32  Y=Al/CLrlt2>  + AL/tULtL)  + A3/CU3tL) 

33  KhlM  "I  HL  AVERAGE  * CK  UM 1 b CK  KhlGhlUEb  1,2, A3  lb” 

3 A L 1 =A  1 * Y/  b 1 + A 1 / U 1 

35  CL=AL*Y/ <bl*bL)  + AL/LL 

36  G3=A3*Y/CbL*b3)  + A3/U3 

37  L ( 1 ) = L 1 / A 1 

3b  V;C2)=C.2/AL 
3 V C 3 ) = £ 3/  A,3 

AO  KhlM  Cl, CL, C 3 

A 1 E=A1*!n(  1 )/A  + AL  + L ( L )/ A + A 2 * '/;  ( 3 ) / A 

A2  KhlM 

A3  KhlM  "IFiE  7C1AL  AVERAGE  V/AI1  l.\  IFiE  bYblEM  lb  < lt\  LAI  b) " 

A A KhlM  to’ 

AS  KhlM 

A 6 A=Al/(UlT3)+AL/(ULt3)+A3/CL3t3) 

A 7 Y=A1/(U1 »L)+AL/ CULtL)+A3/CU3t2) 

Ab  A < 1 ) =2*  A/  b 1 + b*'r*(Al/<bltb)>/(blTb> 

a 9 V ( 1 ) =A ( 1 ) - CWC1)  - ( 1 / U 1 ) ) t 2 + Lit (-2) 

50  AC2)=2*A/ ( (bl tL>*32)  + 2* Y * C Y - A3/ C U3t 2 > > / C C b 1 * 32 ) 1 2 > 

5 1 ACL) =A(2> +L* Y*C A 1/ C U 1 1 2) )/ C C b 1 1 3)*bL> 

5LVCL)=ACL)-CV.CL)-C  l./UL))tL  + 02»C-2) 

53  AC3)=L*A/ C Cb2»2)*b3)  + L* C Y t 2 ) / C C 32* 33 ) t 2 ) 

5 A A(3) =AC  3)+L* Y * C Y-A3/ C U3»  2) )/ < ( bLt  3)*b3) 

55  * 

56VC3)=AC3)- CtoC  3)- C 1 ./l>3)  ) »2  + L3t(-L) 

57  KhlM  "IhE  AV/G  3 Y bl  Eh.  V.All  F L h EAGh  KhlLhllY  lb” 

5 b KhlM  'a(  1 ),L  CL  3) 

5V  KK1M 

6 G KhlM  "IKE  VA  K 1 Ai\  GE  LF  IKE  bYblEK  UlAll  Ft  K EALK  Kbit  hi  YY  lb* 
6 1 KhlM  VC  1 ),  VCD,  VC  3) 

62  KhlM 
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6  3 l MM  "Li  SI  hi  blJ'l  ICt\  b Kh  3 i/1-3  '/.MilhC  IIF-Eb" 

6 4 F L h h = 1 ’1 L 3 

6b  IF  ..(h)  < 40.  CL  '1C  6b 

66  LO\,aO>  = 1 

6 7 CO  1C  179 

6b  IF  V.OO  + CWK)  )t  .b  > AC.  GC  'iL  17  7 

6 9 FLh  1=1  '1C  AG 

7 0 X< 1 ) = 1 . -bC ./ ( ( 10+1 ) t 2) 

7 1 t <h,  1 > = ..  (h)  + ( ( V(h)  t . b)*(  10+1  )/  10.  ) 

72  IF  Y(k,l)  >^G  • CL  10  22 

73  * 

7 b Z(K, I > = 1M ( Y(h,  I ) > + 1 
7 6 b < A 1 - 1 ) = 1 . - X ( I ) 

7 7 1 ( h,  A 1 - 1 ) ( K ) - ( ( V ( h ) t o . b ) * ( 10+  1 )/  1 0 . ) 

lb  IF  1 ( K, I ) > = 0 • GL  10  foO 

79  l(k,l>=0. 

60  U(K,l)  = IM(1<h,I))+l 

bl  NEXT  I 

82  F £ h 1 = 1 10  aG 

b3  i\(k,i)=G 

8 A NEXT  1 

bb  FCh  1 = 1 10  AG 

86  L=U(h,I> 

87  J=A(h,l) 

88  N <h,L)=N(K,L)  + 1 

8 9 t\(k,  J)=i\(K,  J)  + 1 
90  i\EXl  1 

9 1 F.(k,0>=0 

92  h=0 • 

93  F£h  1 = 1 10  AG 
9A  N(K,I>=0 

95  FOE  J=1  10  1 

96  F.(k,  I )=to<K,  I >+N(K,  J> 

97  l\EX1  J 

9 8 h=K(k,  1 ) 

99  L = F <k, 1- 1 ) 

100  FOh  h = 1 10  40 

101  S(A0+N)=X(N) 

102  NEXT  t\ 

103  hCK, 1 )=bCk) -SCL) 

10A  h = F-  + h(k,I) 

10b  NEXT  I 

106  FhlM"  i\ UFifcEh  OF  LAYS"," 

107  0=0. 

106  2 = 0 • 

109  FOh  1=1  10  AO 
1 10  h(K, 1 )=h(K, 1 )/h 
111  G(K,  1 )=F-(h,  1 ) 

1 12  i\EX1  1 
1 13  h=0 

1 1 A FOh  1=1  10  AG 

1 1 b h = h+h(K, 1 ) 
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FhGEAElL 11 Y" 


116  11-  HK*1)>0  GO  1C  1 IV 

1 17  4 = 4 + 1 . 

life  IF  Kh»l+l)>0  GO  1L  12  1 
1 IV  IF  K=1  GC  1C  170 

120  NEX1  1 

121  J=1N1 ( <4+ 1 )/2) 

1 22+ 

12b  IF  Ii\  1 ( 4/  2 ) = 4/  2 GC  10  162 
1 29  + 

130  FGK  1=1  10  40 

131  IF  h(K>  1 >+K(K,  1+  1 )>0.  GO  10  160 

132  * 

1 34  F OK  N = 1 10  J+  1 

13b  G(K.»  I+t\-  1 )=KCK>  I)  / (J+l  ) 

1 36  KEX1  i\ 

1 37  GCK,  1+  1+4)=KK,  1+  1+4)/ ( J + 1 ) 

1 38  FGK  N = 1 1G  J 

140  C(K, I+J+N- 1 )=G(h, l+J+K- 1)+GCK, 1+1+4) 

14  1 + 

14b  NEXT  N 
1 48  + 

1 4V  + 

lbO  GO  10  170 
1 55  + 

160  NEXT  1 

162  FGK  1=1  10  40 

163  IF  K ( K.» I ) + F ( K.* 1 + 1 ) > 0 • GO  10  16V 

164  FOK  N = 1 10  4/2  + 1 

165  G ( I+iSt-1  )=PCKj  I ) / ( J + 1 ) 

166  0(K, 1 + 4/2+N )=FCk, 1 +4+ 1 )/ ( 4/2+ 1 ) 

167  NEXT  N 

168  GG  10  170 

16V  NEXT  I 

170  FOK  1=1  10  40 

171  FKIN1  I- ^'■IHKGUGH'MjGCK,  1 ) 

1 72  G=G+C ( K j I ) 

173  IF  G>l.-10t(-5)  GO  10  175 

174  NEXT  I 

175  FK1M 

176  GO  10  1 80 

177  FOK  1=1  10  40 
176  GCk, 1 )= 1/ 40 
179  KiEXl  1 

160  NEXT  K 

161  blM  C ( 4 1 ) 

165  GC 1 ) = 0. 

190  KK1M 

195  + 

200  FGK  1=1  10  40 

201  0 ( I ) = (Al  + GC I > I ) + A2+G(2,I)  + A3+G<3*l))/A 

202  + 
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2G3  IF  C( I- 1 X 1L. t (-5)  GL  TL  22G 
2Gb* 

2 1G  NEAT  1 
2 1b* 

22G  h=C 

221  FOR  I =1  TG  AO 

222  R=R+ I *C ( 1,1) 

223  NEXT  I 

22a  FRINV’ThE  WE  1 CM  EG  AVERAGE  GF  ThE  PRI-I  GlbTRIbMIGN  lb” 

22b  PRINT  K 

226  PRINT 

227  PRINT 

230  b=  0 

231  GG  10  301 
2 32* 

233* 

234* 

23b* 

236  PRINT 

237  PRINT 

240  T =0 

241  FGR  1=1  TO  40 

242  T =T  + 1 *G ( 3, I) 

243  NEAT  1 

244  PRINT  "THE  WElGhT  EG  AVERAGE  OF  ThE  PRI-I1I  DISTRIBUTION  1SM 

245  PRINT  T 

246  V = 0. 

247  FGR  1=1  TO  40 

24&  v=v+i*c<n 

249  NEAT  I 

250  PRINT 

251  PRINT 

252  PRINT  ”T  RE  WElGhTEG  AVERAGE  FOR  ALL  3 G I 5T  R1  bUT  1GN  b lb” 

253  PRINT  V 
255  GC  TG  3G0 

290  PRINT  ,,HELP-ThE  SYSTEN  IS  G VERLG  AGEG  ! ! ” 

300  bTGP 

301  FGR  I = 1 TG  40 

302  S=S+ I *G( 2,  I ) 

303  NEAT  I 

304  PRINT  ’*  THE  WElGhTEG  AVERAGE  0F  ThE  PR1-II  GI  bT  RI  bUT  I GN  lb” 

305  PRINT  b 

306  PRINT 

307  PRINT 
306  GG  TG  240 
2231* 
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KOLMOGOROV-SMIRNOV  TEST  STATISTIC  FOR  POISSON  FIT 
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1  * Ihlb  15  PROGRAM  PGI5  FCK  lEbllXC  RG155GN  FITS  FG  K 3Aiv  PLE  31AE5 

2 * EEVaEEN  N = 5 i 100  AT  A 5%  LEVEL  GF  5 1 GM  F I CANOE 

3 L1M  P(251)/G<251)/h<251>/5<251)/lC  ICDiF  ( 25 1 ) / C< 25 1 ) / AC  25 1 > 

A N=  10 

5 LATA  35/  lA/A9/53/57/21/30/37/A2/Afc 


6 * 

7 * 

b * 

9 * 

1 0* 

1 5 

Fek  1=1  10 

N 

20 

HEAD  A ( I ) 

25 

NEXT  1 

30 

A = 0 • 

35 

F0R  1=1  10 

N 

AO 

A=A+AC 1 )/N 

A 5 

NEXT  1 

50 

F C 0 ) = 1 • 

55 

FOR  1=1  10 

250 

60 

F C 1 ) = A*F  < I 

- 1 )/  I 

65 

NEXT  I 

70 

P ( 0 ) =EXF  C -i 

A) 

75 

QC0)=P(0) 

80 

FOR  1=1  10 

250 

bS 

P( I ) =P( 0)*F( I ) 

90 

G(1)=G(I-1)+FC1) 

95 

NEXT  1 

100 

k( 1 A)=  1/  10 

1 05 

K(21 ) = 1/  10 

1 10 

RC  30 ) = 1/  10 

1 15 

R<35)  = 1/  10 

120 

h(  37 ) = 1/  10 

125 

hC  A2)  = 1/  1 0 

1 30 

h C 4b ) = 1 / 10 

1 35 

R( 49>  = 1/  10 

1 AO 

k( 53) = 1/  10 

1 A 5 

RC  57 ) = 1/  10 

1 50 

* 

1 55 

* 

160 

* 

165 

* 

170 

* 

1 75 

* 

1 80 

* 

185 

* 

1 90 

* 

1 95 

* 

200 

* 

205 

i>  ( 0 ) =k  ( 0 ) 

210 

Flo  h 1 = 1 10 

250 
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215  SCl)=S(l-n+P<I) 

220  NEX1  1 

225  L(0)=AfcS(C(0)-5C0) ) 

230  FOP  1=1  10  250 
235  b<I)=Afcb(C(I>-S(l>> 

240  NEX1  1 
245  E = 0(0) 

250  FOP  1=1  10  250 
255  IF  D(IXE  GO  10  265 
260  E = L ( 1 ) 

265  NEXT  I 

270  PP1N1  "EPANCh  77^:  C0N1PAC1  MANAGEMENT " 

275  PH  INI  E 

280  LAI  A . 56* • 53* • 50* • 47* • 44* • 4 1 * • 40* • 38* • 37* • 35* • 34*  • 33*  *32*  *31 
285  DATA  . 30*  . 29*  . 29*  . 29*  . 2b*  . 28*  . 27*  . 27*  . 26*  . 26*  . 25*  . 25*  . 24 
290  FOP  1=5  10  30 
295  REAL  1C1) 

300  NEX1  1 

305  IF  N<  = 30  GO  10  330 
310  1 (N>= 1 . 36/ (Nt . 5) 

330  IF  E> 1 (N ) GO  10  350 

340  PPIN1  "YES-lhE  DATA  AhE  POISSON!" 

345  GO  10  400 

350  PRINT  "NO - 1 HE  DATA  APE  NO!  PGISSGN  A1  1 HE  5%  LEVEL!" 

351  DATA  .67* .63*  .60* .56* .53* .49*  *47* .45* .43* .41* .40*  .39* .38 

352  DATA  . 37 * . 36* . 36* • 35* . 34* . 33* . 32* . 32* . 3 1 * • 3 1 * . 30* • 30* . 29 

353  FOP  1=5  10  30 

354  PEAL  1CI) 

355  NEXT  1 

360  IF  N<=30  GO  10  370 
365  1 <N)  = 1 .63/ (Nt . 5) 

370  IF  E > 1 ( N > GO  10  37  5 

372  PPIN1  "bill  IhEY  APE  A1  IhE  1%  LEVEL" 

373  Ge  10  400 

375  PP1N1  "IhE  DATA  APE  STILL  NO!  POISSON!" 

400  STOP 
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1 * THlb  lb  hLuTINE  FLh  TEbTlhG  EL  UAL  1 '1  Y LF  L 1 bT hi bUT 1LN b . 

2* 

3 * K IS  IFt  NUF.bEh  LF  CLAbb  lMEhVALb  Lh  hEbHCNbE  VALULb. 

5 N 1 = 1 7 b 1 
10  N2=N 1 
lb  K= lb 

20  LIN  L ( 30  ) * F.  ( 30  ) , F ( bO  ) , G < 30  ) 

2 1 ' LI F.  1(30) 

2 b LATA  62  3/ a7  b / a£  , a7  / 7 2 / a b / bb  / 6 3/ A b/ b2/ 3 b/ b b / 6 9 / b 1 / b2 
26* 

30  FLh  1=1  TO  k 
3b  hEAL  L ( 1 ) 

AO  F(1)=L(I)/N1 
Ab  NEXT  I 

bO  LATA  2bl, 2bl, 39b, 1AA,  Iaa, 276,  126*69, Al, 27, 23,  12,9 
bl  LATA  7,2 
bb  FLh  1=1  TL  K 

6 0 KEAL  F.  ( I > 

6b  G( 1 )=F < 1 )/N2 
70  NEXT  1 
7b  F ( 0 ) =0 • 
bO  G < 0 ) =0 • 
bb  FCh  1=1  TL  K 
90  F(1)=F(1)+F(I-1) 

9 A* 

9b  G< 1 )=G< I )+G< 1- 1 > 

96  NEXT  1 
97* 

9b* 

1 00* 

10b  FLh  1=1  TL  K 
1 10  L(I)=Abb(F(l)-G(l)? 

1 1 b NEXT  1 

1 1 6 * W E HAVE  JUbT  CEFHUTED  ALL  THE  AbbFLUTE  LIFFEhEb  bET WEN 
1 17*THE  T WE  CUFiULAT  I VE  LlbThlbUTILN  FUNLTILNb. 

1 1 b*WE  FUbT  NEXT  LETEhF.lNE  WHIUh  LF  THEbE  AbbLLUT  E L1F  FEhENCEb 
1 1 9*  1 b THE  LAhGEbT  - WHEN  F,  ULT  IHLI  EL  bY  AN  AHHhL  Hhl  ATE  CLN  bl  AN  T / 

1 20* 1 H 1 b WILL  bE  LUh  TEbT  bTATlbTlt. 

121  E=L<1> 

12b  FLh  1=1  TL  k 

130  IF  L < 1 ) <E  GC  TL  UC 

1 3b  E = L ( 1 ) 

1A0  NEXT  1 

1 Ab  F = <E*<2»<1/2)))/<(1/M+1/N2)t(1/2)> 

IbO  HhlNT  "T HE  VALUE  LF  THE  TEbT  bTATIbTIC  lb"/F 
1 b 1 HrlNT 

1 b2  HhlNT  "THEnEFChE/ " 

1 bb  IF  F> 1 .36  CL  TL  190 

156*1.36  lb  THE  ChlTILAL  b%  hLLF.L  CL  hL  V - bh  1 hNL  V VALUE. 

160  HhlNT  -THE  TWL  F L HuLAT  1 LN  b HAv  E HhLV 1L  EL  hLbHL,\bhb  WHICH  AhE" 

161  Hhll\T  "blAT  IbT  1 CALL  T I bEN  T I CAL  ! " 

16b  bTLF 

190  HhlNT  "T  HE  TWL  L 1 bT  h I HUT  1 LN  b CAN  EE  CLN  bl  LEhF.L  TL  bl" 

191  Hhli\T  "blCislF  ICAnTLY  L IFF  EHFIn  T ! ” 

20C  £lL  F 
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CONVOLUTION  ROUTINE 
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10  * IhlB  FRCChAF.  15  L5LL  10  CGNvGLblE  PROBABILITY  LI51Rlbbll  ON  5 IRRb 
15  * AM  NbFbER  OF  PAlh  ELEFEM5. 

20  * 

25  * THE  FOLLC'/.ING  VAnlABLE5  A*E  b5EL:*‘ 

30  * Nl=NbF.BE'h  OF  bAlh  ELEPEM5 

40  * F 1 1 F2=NbF.BEh  OF  H0bK5  IN  A bAKlICLLAK  LI51R1&U1I0N 

45  * L(I)>E<1)*F(I)  =FhObAbI  L 11  Y L I 51  hi  bUl  ION  AFFECTING  A N ELEFEM 

50  * L C 1 ) = Ii\ bbl  Ll5lKlbbl  ION 

55  * E ( I ) = INI ERNAL  bl  51  nl Bbl  ION 

60  * F ( I ) =0 U1  Pbl  L I 51  hi  bill  I ON 

6 5 * 

70  * 1 HE  FOLLOWING  DATA  LINE5  AhE  KEGblhEL: 

75  * 1000  LAI A Nl 

BO  * 2000  LAI  A hbE(l)#E(2)» * E (F.  1 ) 

90  * KEFEA1  ABOVE  LINE  FGh  ALL  ELEFEN  1 5 

92  * 

9 4 * 

100  DIF.  A(  1C0),D(  100), EC  100), F<  100) 

106  LEI  03=0 
1 10  READ  N 1 
160  REAL  FI 
190  FOR  J=1  10  FI 
200  REAL  L ( J ) 

220  NEXT  J 

230  FOR  1=2  10  Nl 

240  READ  F2 

250  FBR  J=1  10  F2 

260  READ  E ( J ) , 

2 BO  NEXT  J 

290  LEI  F3=F 1 +F2 

300  FOR  J=2  10  F3 

310  LEI  F ( J > =0 

320  FGR  K=  1 - .1  10  FI 

330  IF  J-K< 1 THEN  360 

340  IF  J-K>F2  THEN  360 

350  LEI  F(J)=F< J)+D<K)*E<J-K) 

360  NEXT  K 
362  NEXT  J 

364  IF  I=N1  IhEN  360 

371  F0R  J=I  10  F3 

372  LEI  D(J)=F(J) 

373  NEXT  J 

374  LEI  F1=F3 

375  NEXT  I 
3B0  PRIM 
3B2  PRIN1 
364  PRINT 
390  PRIM 

400  PRIM  "TP.E  NUFBER  OF  bhANChES  IN  1P15  PATH  ARE”;  1 
410  PRINT 
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420  FkIM  ”toAX  LAYS  = ”,*  K3 

4 30  FkIM 

460  H\IM  ’•LAYS”,  "FkG  b-LAY  S" 

470  LSI  Y=  1 - 1 
480  LSI  C 1 =0 
490  FGk  J=1  10  to 3 
300  LEI  G 1 =C 1 +F ( J ) 

3 10  LEI  Y=Y+ 1 

520  IF  Y <8  1 HEN  570 

530  FkIM  J,  F ( J ) 

532  LEI  C2=J/8 

5 36* 

537  LEI  G3=C3+G 1 
539* 

540  LEI  Y = 0 
550  LEI  01=0 
560  GO  10  590 

570  IF  FCJX.001  THEN  590 
580  FkIM  J,F ( J) 

5 9 0~  L E 1 L(J)=F(J) 

600  NEXT  J 

603  IF  Y=8  THEN  610 

605  LEI  C2=G2+ 1 

606  FkIM  '•  ”,  02,  0 1 

608  LEI  G3=  G3+  C 1 

610  FkIM 

6 12  FkIM  " SUto  OF  kkOB-DAY S=”, G3 

7 CO  Z=0 

7 10  FGk  J=I  10  to3 
720  L-L  + J*F  ( J ) 

725  NEXT  J 

726  FkIM 

730  FkIM  ”1HE  AVEkAGE  NUtoEEk  OF  DAYS  IS” 

740  FkIM  Z 
1000  DATA  2 . 

2000* 

2001* 

2002* 

2003* 

2004  DATA  15, . 135, . 174, .075, .092, .09, . 122, .076 

2005  DATA  . 08 1 , . 029, . 028,  . 03,  . 033,  . 0 1 3,  . 0 1 1 

2006  LAI  A .011 
2007* 

2008  DATA  1 4, . 1 59, . 09 1 , . 034, . 053, • 048 , . 1 04, . 1 98 

2009  LAI A . 142, .046, .02, .023, .039, .018, .025 
2010* 

2011* 

2999* 

9999  END 
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EXPONENTIAL  FIT  TEST 
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1 KEF  _lhlS  IS  PktC-KAF.  EXh  Ft  K _1 1 SlJNG_tlX»-t>NLMJ[_AL 

2 KEF  * FIT  S'  FG  h LAN  PL  E'  Vi  AES  bET WEEN" 5_ANL~  SCO  " ‘ 

3 KEP  AT  A 5%  LEVEL  GF  SI  GN 1 F 1 CAN  CE . 

A KLf'.  N WILL  DENOTE  'IhL  TOTAL  NGF.bEh  GF  KIMS. 

5 hEF  LATA  IS  1NPOTEL  FG K T ( 1 ) =”CLGCk”  TIFE0F  ITh  AhhlVAL. 

7 N=  1 0 

6 DIF.  S(501),T(501) 

10  DATA  10,20,3G,40>50,60,70,&0>90> 100 

1 1 KEF 

1 2 KEF. 

1 3 KEF 

1 4 KLF. 

1 5 KLF, 

1 6 KEF. 

1 7 KLF. 

1 8 KEFi 

1 9 KEF 

20  KEF. 

25  FCh  1 = 1 1G  i\ 

30  KEAD  1(1) 

35  NEXT  I 

AO  1(0)=0. 

4 5 F0h  1 = 1 1G  N 

50  S(I)=(N-I+1.)*(T(I)-T(I-1)) 

55  NEXT  1 

60  L= INI (N/ 2 ) 

65  Y 1 =0. 

70  F0K  1=1  1G  L 

75  Y 1 =Y 1 +S( 1 ) 

77  NEXT  1 

80  Y 1 = Y 1/L 

85  Y2=0. 

90  FG  K 1 =L+  1 1G  N 

95  Y2=Y2+S( I ) 

100  NEXT  I 

102  Y2=Y2/(N-L) 

105  G‘= Y 1 / Y2 

109  PhINl  "bhANCh  77  a:  CONThACT  P.ANAGEPlENT  ” 

1 10  PK1N1  G 

115  L1=2*L 

1 16  L2=2* (N-L ) 

120  PK1N1  ”DEGhEES  GF  FKEEDGF  AhE”  LI  "AND"  L2 

121  D 1 =L 1 

122  L2=L2 

125  F2=(D1  + 1 .739)/( . 1 I97*bl  + . 1 108)  - ( D2- 3. 986 ) / ( . 1 A 1 4*0  1 - . 2864 ) 

126  F2=F2+( • 14S-.00  17*0 1 ) + ( . 06 1 5*0 1 - 2 . 706 ) / ( 02+30 • ) 

130  Gl  = (02+1 .739)/( • 1 197*02+.  1 108)  - ( 0 1 - 3 . 986 ) / ( . 1 4 1 A* 02- . 2864) 

131  Gl=Gl+(. 145-. 0017*02)  + ( . 06 1 5*02- 2. 706 ) / ( 0 1 +30 . ) 

135  F 1 = 1 »/G 1 

140  IF  G<F 1 ThEN  200 

145  IF  G> F2  THEN  200 

150  PhlNT  "YES>  IKE  LATA  AKE  EXPONENT  I AL ! " 

175  GO  TG  300 

200  PK1NT  "NG » ThL  LATA  AKE  NOT  EXPONENTIAL!” 

300  ENL 
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